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Preface 


The  purpose  of  this  study  was  to  investigate  the 
effects  of  using  screen  mesh,  flattened  by  rolling,  for  use 
in  screen  regenerative  heat  exchangers.  The  driving  force 
behind  the  concept  is  to  reduce  the  void  volume  of  the 
regenerator.  This  reduced  void  volume  should  improve  the 
overall  performance  of  a  cryocooler.  The  research  was 
sponsored  by  Mr.  Robert  Vacek  of  the  Space  Thermal 
Technologies  Branch  at  Phillips  Laboratory  in  Kirtland  Air 
Force  Base  NM.  The  topic  was  suggested  by  Dr.  Ray  Radebaugh 
at  the  Chemical  Science  and  Technology  Laboratory,  National 
Institute  of  Standards  and  Technology,  Boulder  CO. 

In  the  process  of  preparing  this  thesis  I  had  help  from 
many  people.  In  addition,  to  Mr.  Vacek  and  Dr.  Radebaugh  I 
also  was  blessed  with  the  assistance  of  Mr.  Ron  White  who 
shared  his  years  of  cryocooler  experience  with  me.  Mr.  Tim 
Hancock  and  Mr.  Andy  Pitts  provided  the  technical  assistance 
I  needed  to  get  the  job  done.  Also,  my  advisor.  Major  W. 
Jerry  Bowman,  kept  me  from  straying  from  my  objective  and 
was  always  willing  to  help.  Finally,  this  paper  would  have 
been  unreadable  without  the  persistent  proofreading  of  my 
girlfriend,  Ms.  Kerrie  Lalikos. 

Jeffrey  L.  Wiese 
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Abstract 


This  study  investigated  the  effect  on  heat  transfer  and 
friction  characteristics  for  screen  regenerative  heat 
exchangers  with  the  screen  thickness  reduced  by  rolling. 

The  experiments  were  performed  on  250  and  325  mesh,  304 
stainless  steel  screen  using  helium  gas.  Reynolds  numbers, 
based  on  hydraulic  radius.  Re,  were  between  10  and  100. 

Both  the  Colburn  factor,  5/Pr^^,  and  friction  factor,  /, 
decreased  as  the  screen  thickness  was  reduced.  A 
correlation  for  predicting  /  was  found.  The  coefficient  of 
drag  per  screen,  Cj^,  remained  nearly  unchanged  for 
thicknesses  reduced  not  more  than  30  percent.  The  decrease 
in  Colburn  factor  was  significant  for  Re  less  than  40.  For 
40  <  Re  <  100  the  decrease  in  Colburn  factor  was  less  than 
the  experimental  uncertainty. 
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EXPERIMENTAL  ANALYSIS  OF  HEAT  TRANSFER  CHARACTERISTICS  AND 
PRESSURE  DROP  THROUGH  SCREEN  REGENERATIVE  HEAT  EXCHANGERS 


I .  Introduction 

This  research  experimentally  investigated  the  effects 
of  using  screen  mesh,  flattened  by  rolling,  for  use  in 
regenerative  heat  exchangers.  This  chapter  explains  the 
reasons  for  the  research  and  introduces  the  approach.  The 
background  explains  the  application  of  regenerative  heat 
exchangers  and  how  rolled  screen  has  potential  advantages. 

In  addition,  the  experimental  approach's  history  is 
discussed.  Finally,  the  objective  is  stated  and  the 
approach  used  to  reach  this  objective  is  explained. 

Background 

Regenerators  and  Crvocoolers.  Currently,  Phillips 
Laboratory  at  Kirtland  Air  Force  Base,  New  Mexico  is 
pursuing  using  Stirling  cycle  cryocoolers  for  space  borne 
use.  The  typical  purpose  of  space  borne  cryocoolers  is  to 
cool  infrared  sensors  and  superconductive  electronic  devices 
(Chan,  Tward,  and  Burt  1990il239).  These  applications 
typically  require  regenerative  heat  exchangers. 
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regenerators,  (Chan,  Tward,  and  Burt,  1990:1240).  One 
possible  choice  of  matrix  material  is  wire  screen. 

Walker  provided  definitions  for  czryocoolers  and 
regenerators.  First,  he  defined  a  cryocooler  as  "...  a 
device  or  ensemble  of  equipment  for  producing  refrigeration 
at  temperatures  less  than  120  K"  (Walker,  1983:1).  Later  he 
said  that  regenerators  were  one  of  the  two  typical  types  of 
heat  exchangers  used  in  cryocoolers.  He  stated  "A 
regenerative  heat  exchanger  has  a  single  set  of  flow 
passages  through  which  the  hot  and  cold  fluids  flow 
alternately  and  periodically"  (Walker,  1983:5). 

Walker  provided  a  definition  of  the  ideal  regenerator 
which  also  provided  an  understanding  of  a  regenerator's 
function . 

The  ideal  regenerator  can  be  conceived  as  a 
thermodynamic  "black  box"  accepting  gas  at  temperature 
Tg  and  heating  it  to  Tf,,  After  some  time  the  flow  is 
reversed  and  gas  enters  at  Tf,  leaving  (sic).  The 
pressure  drop  across  the  regenerator  would  be  zero. 

When  used  in  Stirling  engines  where  dead  space  is  an 
important  parameter,  the  ideal  regenerator  has  zero 
void  volume.  (Walker,  1983:36) 

According  to  Walker's  definition,  dead  space  refers  to  the 

volume  filled  by  the  working  fluid  which  is  not  in  the 

compression  or  expansion  space  of  the  cryocooler  (Walker, 

1983:577).  and  T,,  represent  the  cold  and  hot 

temperatures ,  respectively . 

Walker  also  discussed  regenerator  design  for  Stirling 

coolers  as  follows: 


1-2 


for  maximum  heat  capacity  —  a  large,  solid  matrix, 
for  minimum  flow  losses  —  a  small,  highly  porous 
matrix, 

for  minimum  dead  space  —  a  small,  dense  matrix, 
for  maximum  heat  transfer  —  a  large,  finely  divided 
matrix, 

for  minimum  contamination  —  a  matrix  with  no 
obstruction.  (Walker,  1983:44) 

Walker's  third  point,  minimum  dead  space,  is  the 

driving  interest  behind  this  investigation.  The  correlation 

to  dead  space  holds  also  for  non- ideal  cycles.  Tailor  and 

Narayanlchedkar  (1990:1411)  performed  a  nimerical  analysis 

which  allowed  for  imperfect  regeneration  in  a  piston- 

displacer  Stirling  cryocooler.  Their  graphs  illustrated 

that  the  refrigerating  capacity  decreased  as  the  dead  volume 

increased  even  for  non- ideal  regenerators.  In  addition, 

Atrey,  Bapat,  and  Narayankhedkar  (1991:1049)  performed  a 

numerical  analysis  on  Stirling  cycle  wire  mesh  regenerators. 

In  their  analysis  of  mesh  size  (wires  per  inch  of  screen) 

they  varied  the  wire  diameter  which  changed  the  dead  volume 

within  the  regenerator .  They  observed  "...  the  maximum 

effectiveness  is  obtained  at  the  minimum  dead  volume 

combination  for  a  given  mesh  size."  (Atrey,  Bapat,  and 

Narayankhedkar  1991:1049) 

As  the  last  three  references  stated,  decreasing  the 
dead  space  within  a  regenerator  is  an  important  concept.  An 
extension  to  Atrey,  Bapat  and  Narayankhedkar ' s  (1991) 
observation  implies  that  methods  to  decrease  the  dead  space 
for  a  particular  mesh  size  are  worth  considering.  One 


1-3 


possible  method  to  accomplish  this  is  to  flatten  the  screens 
before  assembling  the  regenerator. 

Gary  and  Radebaugh  (1991)  completed  an  improved 
numerical  model  for  calculating  regenerator  performance, 
"Regen  3.1."  The  model  incorporated  six  different  matrix 
geometries.  Wire  screen  was  one  of  the  geometries 
incorporated.  Radebaugh  stated  in  telephone  conversations 
with  the  investigator  (1993)  that  he  would  like  to 
incorporate  rolled  screen  into  his  analysis.  He  had  already 
supezrvised  some  tests  with  200  mesh  screen  that  had  not  been 
fully  analyzed.  In  fact,  he  recommended  the  topic  of 
investing  flattened  screens  for  other  mesh  sizes  to  provide 
breadth  to  the  data  base. 

Thus,  to  approach  the  zero  dead  volume  of  the  ideal 
Stirling  cycle  regenerator,  rolling  the  screen  used  for  a 
regenerator  matrix  provides  possible  advantages.  However, 
pressure  losses  could  prove  to  be  a  problem.  Since  only 
limited  experimentation  with  rolled  screen  regenerators  has 
been  performed  more  data  are  needed  to  gain  an  understanding 
of  its  effect.  Once  this  is  done,  the  feasibility  of  rolled 
screen  in  cryocooler  regenerators  will  be  known. 

Previous  Experimental  Research.  The  research  method 
used  in  this  investigation  is  based  on  three  references. 

The  first  two  references  were  test  results  for  woven  screen. 
The  third  reference  used  unrelated  matrixes,  but  provided 
the  technique  to  evaluate  heat  transfer. 
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The  two  references  for  previous  research  on  woven 
screen  were  Kays  and  London  (1984)  and  Tong  and  London 
(1957).  Kays  and  London  (1984)  compiled  several  series  of 
tests  that  were  performed  in  the  late  1940 's  and  1950 's. 

The  tests  related  to  compact  heat  exchangers.  One  of  the 
matrixes  tested  was  crossed  rod  or  woven  screen.  Tong  and 
London  (1957),  however,  were  the  actual  source  for  the  data 
in  Kays  and  London  (1984).  Tong  and  London  presented 
results  for  hydraulic  Reynolds  numbers  between  5  and 
100,000.  These  two  references  provided  the  comparison  data 
for  this  research.  However,  the  tests  were  performed  under 
different  conditions.  The  working  fluid  was  air  and  the 
geometries  were  much  larger.  These  two  references  also 
provided  the  necessary  theory  to  perform  the  friction 
calculations  for  this  research. 

Pucci,  Howard,  and  Plersal  (1967)  detailed  the  history 
and  application  of  the  single-blow  transient  testing 
technique.  The  matrixes  they  tested  were  not  similar  to 
screen  matrixes.  However  they  explained  how  to  apply  the 
single-blow  transient  testing  technique  to  test  compact  heat 
exchangers.  A  primary  advantage  to  the  single-blow 
transient  testing  technique  is  its  minimal  equipment 
requirements . 

Literature  searches  were  completed  to  find  research 
performed  on  rolled  screen.  No  references  were  found.  The 
searches  were  done  on  the  following  three  data  bases: 
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Defense  Technical  Information  Center  (DTIC),  National 
Aeronautics  and  Space  Administration  (NASA),  and  DIALOG 
Information  Retrieval  Service. 

Objective 

The  objective  of  this  research  was  to  experimentally 
determine  the  advantages  and  disadvantages  of  using  screen 
mesh  flattened  by  rolling  In  regenerative  heat  exchangers. 
There  were  two  aspects  Investigated.  The  first  aspect  was 
the  heat  transfer  properties,  quantified  In  terms  of  number 
of  heat  transfer  units,  and  converted  to  the  Colburn 
factor,  Ji,.  The  heat  transfer  properties  were  expected  to 
Improve  by  rolling  the  screens.  The  second  aspect  was 
friction  losses.  The  friction  loses  were  expected  to 
Increase  with  rolling  the  screen  and  were  quantified  In 
terms  of  the  friction  factor,  /.  The  opposing  effects  of 
heat  transfer  and  friction  losses  provided  the  point  of 
Interest  of  this  research. 

Approach 

The  primary  approach  to  this  research  was  to  build  a 
test  apparatus  that  could  provide  a  step  temperature  change 
by  quickly  changing  the  Inlet  flow  with  the  use  of  a  four¬ 
way  valve.  In  addition,  the  system  measured  Inlet  pressure, 
pressure  drop,  mass  flow  rate,  and  transient  temperature 
response.  The  pressure  and  flow  rate  Information  provided 
the  data  necessary  to  calculate  the  friction,  /,  from 
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equations  In  Tong  and  London  (1957)  and  Kays  and  London 
(1984).  The  transient  temperature  response  was  used  to 
calculate  the  heat  transfer  Information  using  the  single 
blow  transient  testing  technique  (Pucci,  Howard  and 
Plersall,  1967). 


1.  The  test  parameters  were  as  follows: 


1.  Screen  mesh  material:  304  stainless  steel 


2.  Screen  mesh  size:  250  and  325  mesh 


3.  Fraction  roll:  0,  0.15,  0.30,  0.50 


4.  Regenerator  tiiblng  (case)  304  stainless  steel 

5.  Regenerator  dimensions:  15.0  mm  Inside  diameter 

15.9  mm  outside  diameter 
38  mm  long 

6.  Working  fluid:  Helium  gas 

7.  Inlet  pressures  (Pi):  10,  20  and  30  atm 

8.  Reynolds  number  (based  on  hydraulic  radius): 
roughly  10,  20,  40,  and  100 

9.  Initial  temperature:  292  -  297  K 

10.  Final  temperature  (after  step  change):  285  -  300  K 
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II .  Theo 
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rv  and  Experimental  Analvsi 

This  chapter  provides  a  brief  explanation  of  the  theory 
upon  which  the  experimental  technic[ues  to  measure  friction 
and  heat  transfer  were  based.  Then  the  investigator's 
application  of  this  theory  is  explained.  In  addition,  the 
empirical  results  from  previous  research  are  presented. 
Finally,  the  method  of  performing  the  uncertainty  analysis 
is  explained.  Since  a  central  source  of  information  for 
compact  heat  exchangers  is  Compact  Heat  Exchangers.  (Kays 
and  London,  1984)  its  format  and  choice  of  vari£d>le8  were 
followed  as  much  as  possible.  In  the  case  of  the  transient 
temperature  response  some  variables  were  used  from  "The 
Single-Blow  Transient  Testing  Technique  for  Compact  Heat 
Exchangers  Surfaces,"  (Pucci,  Howard  and  Piersall,  1967). 

fyjgtiPh 

Friction  was  analyzed  using  two  different  dimensionless 
quantities.  These  two  quantities  were  the  friction  factor, 
/,  and  the  drag  coefficent,  Cx,.  The  next  few  paragraphs 
describe  both  quantities  and  develop  the  equations  used  to 
calculate  /  and  Q,.  Then  the  empirical  curve  fit  from  prior 
research  is  stated. 

Friction  factor  definition.  Friction  was  one  of  the 
two  major  properties  to  measure.  Kays  and  London  (1984) 
defined  the  mean  friction  factor  (similar  to  Fanning 
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friction  factor  for  flow  through  tubes)  for  a  matrix  surface 


as  /  g  — .  They  applied  conservation  of  momentum  for  a 

heat  exchanger  with  uniform  frontal  and  exit  areas  to  obtain 
the  following  equation. 


V, 


'^gcPx 


+f-r— 

A  V, 


(2.1) 


where 

P  = 

= 

G  = 
Sc  “ 

w  = 

AP  = 


density  (inverse  of  specific  volume  v) 

unit  surface  shear  stress 

exchanger  flow-stream  mass  velocity  {W/A^) 

proportionality  factor  in  Newton's  second  law 

mass  flow  rate 

pressure  drop  across  matrix 

pressure  at  entrance  to  the  heat  exchanger 

specific  volume  at  exchanger  entrance 

specific  volume  at  exchanger  exit 

mean  specific  volume 

porosity  of  matrix  surface 

exchanger  total  heat  transfer  area 

exchanger  minimiun  free-flow  area  {pA^.) 

exchanger  total  frontal  area 


Note:  The  first  term  inside  the  scpiare  brackets  is 
flow  acceleration  and  the  second  term  is  the  core  friction. 


Friction  Operating  Equation.  The  investigator  solved 
Eq.  2.1  for  the  friction  factor. 


V,  A 


'G\ 


-(• 


(2.2) 


The  next  step  was  to  apply  simplifying  assumptions  to 
make  this  equation  a  function  of  the  various  measured 
parameters.  The  first  assumption  was  that  the  gas  was 
thermally  perfect.  Next,  the  mean  specific  voliime  was 
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approximated  as  =  (vj  +  v2)/2.  This  approximation  for 
was  recommended  by  Kays  and  London  (1984).  (The  specific 
method  used  to  calculate  the  screen  properties  (Af  A^,  p, 
and  r,,)  and  mass  flow  rate,  W,  are  explained  in  later 
sections  of  this  chapter) .  In  addition,  frontal  area,  A^., 
was  assumed  circular.  The  measurements  were  made  in  metric 
units  so  =  1.  In  addition,  the  effects  of  Mach  number 
were  assumed  negligible  since  the  Mach  numbers  were  less 
than  0.1.  The  testing  was  near  ambient  temperature  and  the 
regenerator  case  was  insulated.  Thus,  the  process  was 
assumed  adieOaatic  for  the  friction  measurement  experiments. 
These  last  assumptions  allowed  the  static  and  stagnation 
temperatures  to  be  the  same  and  Tj  »  Thus,  the  final 

operating  equation  for  the  friction  factor  was 


D  was  the  inside  diameter  of  the  regenerator  tube  and  R  was 
the  gas  constant  for  helium. 

Previous  Results.  Kays  and  London  (1984)  provided  a 
graph  of  /  vs.  Re  for  an  infinite  randomly  stacked  woven- 
screen  matrix.  This  graph  was  an  excerpt  from  Tong  and 
London  (1957).  Tong  and  London  provided  curves  of  /  vs.  Re 
for  lines  of  constant  porosity.  In  addition,  they  presented 
graphs  of  Cp  vs.  Re^  that  were  the  basis  of  their  empirical 
curve  fit. 


Re 


-0.33 


0.54 

P 


(2.4) 


where 

Cq  =  drag  coefficient  per  screen 

Re  =  Reynolds  niimber  for  tube  like  flow  or  internal 
flow 

Rcj  =  Reynolds  number  for  flow  over  a  body  or  external 
flow 

Note:  Both  Reynolds  niunbers.  Re  and  Re^,  are  explained 
in  further  detail  in  a  later  section  in  this  chapter. 

Tong  and  London  derived  the  following  conversion 
relationship  for  "essentially  constant-density  flow"  (Tong 
and  London,  1957:1561): 


(2.5) 


where 

6  ==■  screen  thickness 

=  hydraulic  radius,  LAjA 

L  =  Length  of  matrix  in  flow  direction  (tube  length  - 
gap  (explained  in  Chapter  III)) 
a  -  ratio  of  minimum  free-flow  area  to  frontal  area 

a  =  {X-d-meshY 
d  =  wire  diameter 
mesh  =  mesh  size,  wires/inch 

Eq[uation  2.5  proved  very  useful  when  converting  test  data 
for  /  to  Cp  which  allowed  another  perspective  as  well  as  a 
comparison  to  the  results  of  Tong  and  London. 


Heat  Transfer 

Common  practice  calls  for  the  heat  transfer  results  to 
be  in  the  form  StPr^.  This  term  will  be  commonly  referred 
to  in  its  shorter  form  as  the  Colburn  factor,  j„=StPr^'’. 
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First,  the  method  detailed  by  Pucci,  Howard  and  Piersall 
(1967)  to  find  the  dimensionless  quantity  number  of  heat 
transfer  units,  N^,  is  explained.  Next  this  section 
explains  the  transition  of  to  Finally,  the 
investigator's  application  of  the  entire  process  is 
detailed. 

Sinole-Blow  Transient  Testing  Technique.  The  single 
blow  transient  testing  technique  used  in  the  investigation 
was  based  on  the  writing  of  Pucci,  Howard  and  Piersall 
(1967).  They  summarized  the  historical  development  of  the 
technique  by  other  researchers.  They  also  added  their  own 
modification.  The  following  is  a  direct  quotation  from 
their  paper  explaining  the  final  result: 

....  The  analysis  is  based  upon  an  energy  balance 

on  an  element  of  the  porous  solid  . . . 

Assiimptions  made  in  the  analysis  are: 

(a)  Properties  of  the  fluid  are  temperature 
independent 

(b)  Fluid  flow  is  steady 

(c)  Porous  solid  is  homogeneous 

(d)  Thermal  conductivity  of  both  fluid  and 
solid  is  infinite  perpendicular  to  the 
fluid  flow  direction 

(e)  Thermal  conductivity  of  fluid  is  zero  in 
the  flow  direction. 

Initial  boundary  conditions  are: 

(a)  The  matrix  is  initially  at  a  uniform 
temperature 

(b)  At  time  equal  to  zero,  the  temperature 
of  the  entering  fluid  changes 
instantaneously  to  a  different,  constant 
value,  i.e.,  a  step  change  in  fluid 
temperature 

(c)  The  matrix  boundaries  are  adiabatic 

(Pucci,  Howard  and  Piersall,  1967:29) 
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One  additional  assumption  was  that  the  thermal  capacity 
of  the  gas  (mass  x  specific  heat)  contained  in  the  matrix 
was  much  less  than  the  thermal  capacity  of  the  solid. 

Thermal  conduction  in  the  solid,  parallel  to  the 
direction  of  the  flow  was  represented  by  the  dimensionless 

k  A 

longitudinal  conduction  parameter  defined  as  A  =  ~  ^  . 
where 

kj  =  matrix  thermal  conductivity 

=  solid  matrix  cross-sectional  area  available  for 
thermal  conduction 
L  =  length  of  matrix 

Cp  =  specific  heat  at  constant  pressure 
Approximation  of  Zero  Longitudinal  Conduction.  Since  A 
was  an  important  factor  in  the  procedure,  an  estimate  of  A 
was  required  before  proceeding  with  further  calculations.  A 
simple  and  straight  forward  approximation  of  thermal 
conductivity  perpendicular  to  the  screens  for  liquid 
saturated  wrapped  screen  is  (Chi,  1976:50) 

where 

k^  =  effective  thermal  conductivity 
ki  =  liquid  thermal  conductivity 
p  =  porosity 

The  approximation  for  A  resulted  by  substituting  the 
thermal  conductivity  of  helium  for  and  realizing  that 
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m  For  constant  length  and  frontal  area  the  maximiun 

value  for  occurred  with  minimum  mass  flow  and  maximum  k, 
(minimum  porosity).  As  shown  in  Appendix  C,  p^„  =  0.4245 
which  yielded  k^  =  0.524  W/mK.  In  addition,  the  frontal 
area  was  0.176  x  10~^  m^,  the  length  was  0.0371  m  (Appendix 
C),  and  the  minimum  flow  rate  was  approximately  0.28  x  10”3 
kg/s  (Appendix  D) .  Thus,  *  0.0017.  Pucci,  Howard  and 
Piersall  (1967:30)  tabulated  finite  difference  results  that 
included  A  but  the  minimum  non-zero  value  was  0.005  which  is 
greater  than  the  estimated  value  for  A^.  Therefore,  the 
investigator  used  the  approximation  of  zero  thermal 
conductivity  in  the  direction  of  the  flow. 

The  approximation  A  «  0  allowed  the  use  of  the  exact 
solution  below. 


slope  = 


(2.7) 


where 


M 

0  = 
-U,(ix)  = 


initial  temperature  of  matrix  and  fluid 

inlet  temperature  of  fluid  after  time  equal  zero 

outlet  temperature  of  fluid 

time  parameter,  hAQ/(W/!g) 

number  of  heat  transfer  units,  hA/ {Wcp) 

convection  heat  transfer  coefficient 

mass  of  matrix 

matrix  material  specific  heat 
time 

modified  (Hyperbolic)  Bessel  Function 
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Equation  2.7  predicts  slope  of  the  dimensionless 
temperature  as  a  function  of  and  time.  From  the  equation 
the  slope  can  be  found  for  any  time. 

As  explained  in  Pucci  Howard  and  Piersall  (1967);  for  a 
given  the  maximum  slope  (using  Eq.  2.7)  has  a  unique  value. 
Thus,  if  the  maximum  slope  can  be  measured  experimentally  the 
can  be  found. 

The  maximum  slope  for  a  given  was  found  by  gradually 
increasing  the  time  parameter,  t.  The  slope  steadily 
increased  until  the  maximim  value  was  reached.  Once  the 
slope  began  to  decrease,  the  maximum  slope  for  the  given  was 
known.  Pucci,  Howard  and  Piersall  performed  this  proceedure 
for  up  to  60.  This  was  insufficient,  so  the  investigator 
extended  the  curve.  The  resulting  vs.  maximum  slope  is 
graphically  illustrated  in  Fig.  2.1  and  tabulated  in 
Appendix  B. 
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Application  of  Single  Blow  Technique.  The  investigator 


defined 


r*  = 


In 


(2.8) 


In  addition,  the  dimensional  values  which  defined  and  r 
were  svibstituted  into  the  left  hand  side  of  Eg.  2.7  to 
obtain 


WCj,  d0 


(2.9) 


For  each  test,  the  maximum  slope  defined  by  Eg.  2.9  was 
found  from  the  experimental  values  of  temperature,  time, 
screen  mass,  mass  flow  rate  and  tabulated  specific  heats. 
Knowing  the  maximum  slope ,  the  could  be  determined  from  the 
vs .  maximum  slope  graph  in  Fig  2 . 1  or  the  ted}le  in  Appendix 
B.  The  next  step  was  to  consider  the  definition  of  Stanton 


number,  St=— — .  From  the  definitions  of  G  (Eg.  2.1)  and 

GCp 

Stanton  number,  5!r,  the  desired  Colburn  factor,  jfj,  can  be 
related  to  N^. 


(2.10) 


Method  to  Find  Maximum  Slope .  The  method  the  investigator 
used  to  find  maximum  slope  relied  heavily  upon  the  spreadsheet 
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Quattro  Pro  SE.  As  described  in  Chapter  III,  the  transient 
temperature  response  was  recorded  to  data  files.  For  each 
test  the  investigator  would  read  the  data  files  into  a  base 
spreadsheet  that  was  already  configured  for  the  analysis. 

dr* 

The  goal  was  to  find  — —  as  required  from  Eq.  2.9.  To 

max 

achieve  this  goal,  the  investigator  made  a  slight 
modification  to  the  equation  for  T*  due  to  the  slight 
difference  (less  than  0.3  K)  that  sometimes  existed  between 
measured  values  for  Tj  and  steady-state  before  the  step 

temperature  change.  The  following  approximation  was  used 


where 


(2.11) 


Tf\  -  average  inlet  temperature  after  the  step  change 

Ti\  =  average  inlet  temperature  before  the  step  change 

Tn  =  average  outlet  temperature  before  the  outlet 
temperature  began  to  change. 

A  quick  summary  using  Fig.  2.2  is  useful.  The  two 
curves  and  7^  represent  the  measured  temperatures  at  the 
inlet  and  outlet,  respectively,  of  the  regenerator.  The 
horizontal  distances  shown  with  r/i,  Jn,  and  Tn  illustrate 


the  time  interval  for  which  the  average  values  were 
calculated.  Equation  2.11  was  applied  to  Tf\,  Tn,  and  Tn  to 


calculate  T*„u!asured'  Finally,  the  heavy  line  labeled  T*re^ssion 
(scale  is  on  right  hand  side  of  graph)  was  the  result  of  a 
linear  regression  performed  on  rVeowred  the  time  interval 
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Time  (seconds) 


dT* 

Figure  2.2.  Sample  Quattro  Pro  Analysis  of  -  Showing 

dO  „„ 

Tf\ ,  Tn ,  and  Tn 


Before  proceeding,  an  explanation  Tf\  is  necessary. 
Since  the  actual  inlet  temperature  was  not  a  true  step 
function,  the  investigator  made  an  approximation.  The 
approximation  was  based  on  an  average  temperature.  T/i  was 
the  average  of  Tj  over  the  time  interval  beginning  when  T, 


had  changed  roughly  1  K  and  ending  at  the  time  of 

(maximum  slope  of  r^) . 

Once  the  values  for  Tf\,  Tn,  and  Tn  were  found  the 
process  of  finding  the  slope  began.  First,  these  values 
were  substituted  into  Eq.  2.11  to  plot  T*.  From  the  plot  of 
T*  vs.  time,  a  time  interval  was  visually  selected  to 
perform  a  linear  regression  of  T*  vs.  time.  Then  the  linear 
regression  result  was  plotted  with  T*  to  ensure  a 
satisfactory  regression  had  been  selected.  If  the 
regression  appeared  to  be  unsatisfactory  the  time  interval 
was  altered  appropriately  and  its  plot  inspected.  This 
process  was  continued  (using  the  same  T'/i,  Tn,  and  Tn)  until 

the  linear  regression  accurately  represented 


Finally,  the  maximum  slope  was  calculated  using  Eq.  2.9. 
Then  the  investigator  used  the  curve  fit  from  Appendix  B  to 


calculate 

=  a(maximum  slope)^ 


maximum  slope _ a _ b 

[0.92,  2.00)  11.77  2.084 


[2.00,  2.80)  12.19  2.023 

[2.80,  5.32)  12.38  2.008 


(2.12) 


Equation  2 . 12  allowed  automation  of  the  data  analysis 
in  the  same  spreadsheet  used  to  calculate  the  friction 
results . 

Finally,  the  Colburn  factor,  jfj,  was  calculated  using 
Eq.  2.10. 
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Previous  Results.  Pucci,  Howard,  and  Piersall  (1967) 
performed  their  tests  on  matrixes  other  than  woven  screen. 
Hovrever,  Kays  and  London  (1984)  provided  a  graph  of  StPr^^ 
vs.  Re  for  infinite  randomly  stacked  woven-screen  matrix. 
This  graph  was  an  excerpt  from  Tong  and  London  (1957).  Tong 
and  London  provided  curves  of  StPr^^  vs.  Re  for  lines  of 
constant  porosity.  In  addition,  they  provided  an  empirical 
cuzrve  fit  for  screen  or  crossed  rod  matrixes  and 
0.55  <  p  <  0.85. 

J„  =  StPr^'*  =  0.375  Re*^  (2.13) 

The  two  equations  below  were  iterated  to  calculate  the 
modified  Reynolds  number.  Re'.  The  first  iteration  started 
by  letting  Re'  =  Re.  Then  the  first  iteration  for  F  was 
found  and  the  second  iteration  for  Re'  was  calculated.  The 
solution  converges  after  several  iterations. 

F=  1.155 - 0.060 l(log,o Re'),  Re’=^-^Re  for  Re'  <1800 

Fp 

Screen  Properties 

This  section  details  the  assumptions  and  equations  used 
to  calculate  the  screen  properties  (A,  p,  A^,  r^,  and  R) . 

Matrix  Total  Area.  A.  The  primary  assumption  by  the 
investigator  was  that  the  area.  A,  was  approximately  equal 
to  the  area  of  a  long  thin  wire  with  diameter  d  and  mass  W,. 
In  addition,  effects  of  the  regenerator  case,  or  tube,  were 
assumed  negligible,  since  the  tube's  surface  area  was  small 
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compared  to  the  surface  area  of  all  the  screen  it  contained. 
Thus,  the  resulting  equation  was 


where 


(2.14) 


-  measured  mass  of  the  screen  in  the  regenerator 
d  -  wire  diameter  as  published  by  the  manufacturer 
=  density  of  solid  304  stainless  steel  (Appendix  A) 

Porosity.  The  porosity,  p,  of  the  matrix  was 

calculated  using  masses.  The  investigator  derived  Eq.  2.15. 

He  assumed  the  matrix  to  be  homogeneous  and  used  the  ratio 

of  the  measured  mass  to  the  mass  of  a  solid  cylinder. 


(2.15) 


/!D^P.L 

Exchana«i»7r  Minimum  Free  Flow  Area.  The  exchanger 
minimum  free  flow  area,  A^,  was  derived  from  an  equation  in 
Kays  and  London  (1984:44), 


(2.16) 


Hydraulic  Radius.  The  hydraulic  radius,  r^,  was 
calculated  from  another  relationship  from  Kays  and  London 
(1984:44) . 

A 


(2.17) 


Fraction  Roll.  The  investigator  derived  a  relationship 
to  calculate  the  fraction  roll,  R. 


(2.18) 
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where 


6  -  measured  thickness  of  screen 

-  measured  thickness  of  screen  before  rolling 

Mass  Flow  Rate 

As  will  be  discussed  in  Chapter  III/  the  instrument 
used  to  measure  mass  flow  rate  was  a  rotameter  calibrated 
for  standard  cubic  feet  per  minute  (SCPM)  of  air.  Since  the 
rotameter  was  calibrated  for  air  the  investigator  needed  a 
method  to  convert  the  readings  to  heliim.  Dally,  Riley  and 
McConnel  (1984)  provided  a  description  of  the  theory  of 
rotameters.  A  rotameter  is  a  tapered  tube  with  a  "bob" 
which  is  suspended  in  the  flow  of  the  fluid  to  be  measured. 
Because  of  the  shape  of  the  "bob",  the  effect  of  viscosity 
was  negligible.  Dally,  Riley  and  McConnel  (1984)  derived 
the  following  linear  relationship,  because  of  the  geometry 
of  the  tube  and  bob. 

^=K^(pb-Pf)pfy  (2.19) 

where 

=  mass  flow  rate  of  the  flow 
K  =  constant  based  on  geometry 
Pi,  -  density  of  the  bob  (metal) 

Pf  =  density  of  the  flow 
y  =  height  of  "bob"  in  tube 

The  investigator  used  Eq.  2.19  to  derive  the  conversion 
to  mass  flow  rate  of  helium.  Since  and  p,,gu„„,  are  both 
much  less  than  pj,  the  ratio  of  flow  rates  for  two  flows 
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with  Identical  rotameters  with  the  bob  in  the  same  position 
is 


_  IPh* 
Kr  iPa^  * 

where 


(2.20) 


Wfj^  =mass  flow  rate  of  helium  at  test  conditions 
=  density  of  helium  at  test  conditions 
Wg^  =mass  flow  rate  of  air  at  standard  conditions 
=  density  of  air  at  standard  atmosphere 

The  mass  flow  rate  of  the  air  can  be  found  from: 

=  (2.21) 

The  calibration  plate  on  the  rotameter  stated  that 

SCFMgi^  =  0.1345  X  percentage  scale  of  the  meter. 

Furthermore,  the  plate  defined  a  standard  atmosphere  as 

70  and  1  atm.  Using  this  relation  and  assuming  thermally 

perfect  gases  the  investigator  derived  the  equation  to 

calculate  helium  mass  flow  rate  from  the  rotameter  reading. 

W=W„^  =  1.526  X  10“®  X  (%scale)  kg/  s  (2.22) 

where 


^iiaro~  barometric  pressure  in  Pa 

P _ gage  pressure  of  the  flow  entering  the  meter  in  Pa 

Ti  =  temperature  entering  the  regenerator  in  °K 

The  temperature  was  assumed  to  remain  constant  since 

all  tubing  essentially  remained  at  room  temperature.  This 

was  accepted>le  because  all  flow  rates  were  measured  before 

the  step  temperature  change  was  initiated. 
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Reynolds  Ninnber 

There  were  two  types  of  Reynolds  number  used  in  this 
investigation.  The  internal  flow  Reynolds  number.  Re,  was 
used  to  present  data  in  the  standard  format  used  in  Kays  and 
London  (1984).  At  times  however,  a  Reynolds  number  based  on 
external  flow,  Re^,  was  found  to  be  useful  in  the  discussion 
in  Chapter  IV. 

The  first  Reynolds  number  was  for  internal  flow.  Kays 
and  London  (1984)  presented  the  definition 

4nG 


(2.23) 


where 


fV  w 

G  =  flow  stream  mass  velocity  ( —  or  - ) 

4  P^fr 


=  hydraulic  radius  (— i) 

A 

H  -  viscosity 

The  viscosity  is  listed  in  Appendix  A  and  was  based  on 
the  regenerator  inlet  temperature  T;. 

The  second  type  of  Reynolds  number  was  for  external 
flow.  Tong  and  London  (1957)  presented  the  definition 


RCd  = 


(2.24) 


where 


d  wire  dieuneter 

maximum  mass  velocity  in  the  matrix  (W/aA^) 
a  =  ratio  of  minimum  free-flow  area  to  frontal  area 
a  =•((*,-  l)/{x,)}^  =  (1  -  dmesh)^ 

X,  =  ratio  of  transverse  spacing  to  wire  diameter 
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Uncertainty  Analysis 

An  uncertainty  analysis  was  performed  on  mass  flow 
rate,  Reynolds  niunber(s),  friction,  and  heat  transfer. 
Different  approaches  were  used  for  each  parameter,  as 
appliced^le. 

Mass  Flow  Rate.  As  explained  earlier  in  this  chapter, 
the  mass  flow  rate  was  calculated  using  Eq.  2.20.  The  form 

fp  Tp 

of  the  equation  was  W = C^-'^  ^  x  (%scale) ,  where  C  was  a 

constant.  The  inyestigator  found  that  ^ meter >  ^7  were 

accurate  compared  to  %scale.  Due  to  the  limited  resolution 
of  the  scale  on  the  rotameter,  the  inyestigator  assiimed  the 
resolution  to  which  the  scale  could  be  read  was  0.5  %scale 
(i.e.  i!t.{%scale)  -  0.5  %).  So,  the  method  of  uncertainty 
analysis  presented  by  Beckwith,  Buck,  and  Marangoni 
(1982:269)  was  applied  with  %scale  as  the  only  yariable.  The 
uncertainty  of  mass  flow,  W,  is  dtW. 


AW  _  I  ^  A(yoscale)  dW  V  _  A{%scale)  _ 

W  W  t5(%scale))  %scale  '  *  ^ 

Reynolds  Number.  The  two  Reynolds  numbers  were  defined 
in  Eqs.  2.23  and  2.24.  The  properties  of  the  matrix  {A, 

Ap,  L,  a,  and  p)  appeared  to  be  accurate  in  comparison  to 
the  resolution  of  the  rotameter.  So,  using  the  same 
approach  as  shown  for  mass  flow  rate  the  uncertainty  of  Re 


ARe  ARcj  A{%jco/^) 

~  —  •  I  A  •  ^  O  j 

Re  Re^  %scale 

Friction.  As  shown  earlier  in  this  chapter,  the 
friction  factor,  /,  was  calculated  using  Eq.  2.3.  The  two 
variables  that  were  of  most  concern  were  the  differential 
pressure,  AP  and  the  mass  flow  rate,  W.  The  uncertainty 
was  found  by  using  a  spreadsheet  to  vary  AP  by  its 
uncertainty,  A(AP) ,  and  vary  W  by  AW.  The  experience  gained 
while  calibrating  the  differential  pressure  transducer 
provided  the  reference  for  its  accuracy.  To  find  the 
uncertainty  in  /,  the  voltage  used  for  calculating  AP  was 
varied  by  ±0.0005  volts  and  the  rotameter  reading  was 
varied  by  ±0.5  %scale.  The  conclusion  was  that  the 
uncertainty  in  /  ranged  from  12  percent  at  the  lowest  Re  and 
decreased  to  1.1  percent  at  the  highest  Re . 

Heat  Transfer.  The  heat  transfer  uncertainty  was 
calculated  using  the  method  of  Beckwith,  Buck,  and  Marangoni 
(1982).  A  rigorous  solution  would  involve  partial 
derivatives  of  modified  Bessel  functions.  Instead  the  curve 
fit  derived  earlier  was  employed  since  it  is  accurate  in  the 
region  of  interest . 

As  explained  earlier,  the  heat  transfer  was  calculated 


using  Eq. 


2.10  (J„=St  )  • 


Again,  assuming  the 


matrix  properties  were  accurate  and  applying  the  curve  fit 
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=  a{max slope)’’  to  Eq,  2.9,  the  resulting  ecjuation  was  found 
to  represent  the  uncertainty  in 


=  b 


Jh 


(dT* 

\ 

A 

[  dO 

mix  > 

dT* 

de 

mix 

(^W\ 


(2.27) 


The  next  step  was  to  approximate  A 


(dT* 

\ 

mix  > 

This  was 


accomplished  by  utilizing  the  statistical  output  which 
Quattro  Pro  SE  generated  with  the  linear  regression  used  to 
dT*\ 


calculate 


de 


The  percentage  error  of  the  slope  was 


calculated  by  dividing  the  "standard  error  of  the  x- 
coefficient"  by  the  "x-coefficient"  (Borland  1991:479). 
However,  performing  the  ratio  for  each  test  was  not  as 
revealing  as  evaluating  this  ratio  for  a  random  ssunple  and 
finding  the  upper  bounds  of  the  percentage  error.  In 
addition,  the  investigator  observed  that  the  percentage 
error  was  dependent  on  the  mass  flow  rate  of  the  helium.  As 
a  result,  the  percentage  error  was  calculated  for  a  random 
sample  of  five  tests  for  each  of  the  four  different  mass 
flow  rates  used  in  the  experimental  tests.  The  result  was 
that  the  bounds  were  as  follows: 


W  (10”3  kg/s) 

0.28 

0.56 

1.2 

2.9 

bound ( %error) 

4 

8 

13 

25 

A  linear  regression  of  the  bounds  resulted  in  the  equation 
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7,700Fr+3.07%. 


Thus,  the  final  equation  for  estimating  the  Colburn 


was 


^  =  3^(77  fF+ 0.0307 f 


(2.28) 

factor 


(2.29) 


2-21 


Ill .  Experimental  Design 


In  this  chapter  the  equipment  used  for  the  experiment 
is  described.  First,  the  overall  test  apparatus  is 
explained  and  illustrated.  Then,  part -by-part  each 
component  of  the  apparatus  is  explained.  This  will  help  to 
explain  the  reasons  behind  any  decisions  and  provide  a 
detailed  description  of  each  component  and  method  of  use. 

The  effect  of  components  on  the  data  and  results  is 
discussed  in  Chapter  IV. 

Overall  Design 

As  explained  in  Chapter  II  an  apparatus  was  required 
that  provided  pressure  and  temperature  sampling  before  and 
after  the  flow  entered  the  regenerator.  In  addition,  it 
needed  to  produce  a  step  temperature  input .  Figure  3.1  is  a 
schematic  of  the  entire  apparatus,  and  Fig.  3.2  is  a 
photograph  of  the  entire  apparatus. 

Working  from  left  to  right  across  Fig.  3.1  the  first 
item  is  the  two  helium  tanks  which  provided  the  helium  for 
the  system.  Then,  a  pressure  regulator  controlled  the 
pressure  of  the  system.  After  the  pressure  regulator  the 
flow  was  split  to  provide  a  cold  and  warm  flow  path.  The 
warm  flow  was  room  temperature  and  the  cold  flow  was  chilled 
with  an  ice  bath.  Both  flows  nearly  crossed  paths  at  the 
four-way  valve  which  was  quickly  switched  when  the  step 
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Figure  3.1.  Test  Apparatus  Schematic 


temperature  change  was  required  (shown  in  the  initial 
position) .  Then  the  flows  split  into  the  primary  and 
secondary  flows.  The  primary  flow  actually  flowed  through 
the  regenerator  for  the  tests.  The  inlet  and  exit 
temperatures,  Tj  and  T2>  were  measured  as  well  as  the  inlet 
pressure,  Pj,  and  the  pressure  drop,  aP.  The  primary 
metering  valve  provided  precise  mass  flow  control  through 
the  regenerator.  Finally,  the  flow  meter  measured  the  flow 
rate  through  the  regenerator.  The  secondary  flow  allowed 
the  cold  flow  to  chill  the  lines  leading  to  the  four-way 
valve.  The  temperature  leaving  the  four-way  valve  was 
measured  through  the  secondary  flow  line.  A  metering  valve 
controlled  the  flow  rate  through  the  secondary  flow  tubing 
(kept  to  a  minimum) .  The  two  flows  recombined  in  the  vent 
line  which  vented  to  an  exhaust  fan. 

Regenerator 

The  regenerator  design  was  based  on  a  design  used  by 
Rawlins  (1992)  at  the  Chemical  Science  and  Technology 
Laboratory,  National  Institute  of  Standards  and  Technology 
(NIST) .  These  tests  were  performed  with  the  guidance  of 
Radebaugh  who  recommended  investigating  rolled  screen  as 
mentioned  in  Chapter  I.  The  NIST  regenerator  tube  was  made 
of  stainless  steel  with  a  15.875  mm  (5/8  inch)  outside 
diameter,  a  0.254  mm  (0.010  inch)  wall  thickness  and  a 
length  of  122  mm  (8.80  inch).  Figure  3.3  shows  the  design 
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for  this  thesis.  The  regenerator  was  the  same  outside 
diameter  but  the  length  was  significantly  less  than  the  NIST 
regenerator.  The  shorter  length  reduced  the  quantity  of 
screen  and  labor  required.  In  addition,  the  lower  pressure 
drop,  AP,  could  be  measured  by  pressure  transducers  "on- 
hand”  at  the  AFIT  Department  of  Aeronautics  and  Astronautics 
Laboratories . 


Figure  3.2.  Photograph  of  Test  Apparatus 


As  mentioned  in  Chapter  I,  250  and  325  mesh  screen  was 
used  so  that  the  results  could  augment  earlier  results  at 
NIST.  The  NIST  tests  were  performed  on  200  mesh  screen  (the 
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NIST  results  were  not  analyzed  in  time  to  be  included  in 
this  research).  The  actual  material,  304  stainless  steel 
was  selected  because  of  screen  availability.  Both  screen 
mesh  sizes  were  Buffalo  Wire  Cloth  manufactured  by  Buffalo 
Wirework  Company,  Inc.  The  250  mesh  was  plain  weave  and  the 
325  mesh  was  twill  weave. 


5/8  in.  (15.879  mm)  304  Stainless  Steel  Tube 
Approximately  1.5  in.  (38  mm)  long 


Thermocouple 

Connection 


2  Oversized  Screens 
5/8  in.  Diameter 
Edoxv  (15.879  mm) 


Flow 


Location  of 
Thermocouple 
Junction 


Lead  Solder 

u  u 

Pressure  Port 

1/8  in.  (3.175  mm)  Copper  Tubing 


Modified  Swagelok  Reducing  Union 
5/8  in.  to  3/8  in. 

(15.879  mm  to  9.525  mm) 

Figure  3.3.  Regenerator  Schematic 


As  shown  in  Fig.  3.3,  the  regenerator  design  utilized 
"off-the-shelf"  hardware.  The  design  provided  a  simple  yet 
easily  removed  and  replaced  regenerator.  Swagelok  brand 
fittings  provided  strength  and  easy  installation.  In 
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addition,  the  regenerator  housed  the  temperature  and 
pressure  sample  points. 

Screen  Preparation.  Preparing  the  screen  was  the  most 
difficult  and  time  consuming  task  in  making  the  experimental 
apparatus.  The  difficulty  was  primarily  due  to  the 
unavailability  of  equipment  on  Wright  Patterson  Air  Force 
Base  OH,  designed  specifically  to  roll  material,  such  as 
screen,  to  the  precision  and  repeatability  required  for  this 
research.  The  time  requirement  resulted  from  the  tedious 
tasks  of  rolling,  cutting,  cleaning  and  packing  nearly  5000 
screen  disks.  The  next  few  paragraphs  describe  the  method 
and  equipment  used  for  these  tasks. 

Screen  Rolling.  As  mentioned  above,  no  equipment 
was  readily  available  to  roll  the  screen  in  a  precise  and 
repeatable  manner.  However,  after  much  experimentation,  a 
Milwaukee  Horizontal  and  Vertical  Milling  Machine,  model  205 
S-12,  proved  satisfactoi^.  The  machine  was  used  in  the 
horizontal  mode. 

The  moving  table  of  the  mill  was  the  key  to  this 
application.  A  standard  1  inch  x  4  inch  x  24  inch  hardened 
machinist  parallel  bar  was  attached  to  the  moving  table  and 
used  as  the  rolling  surface  on  which  the  screen  was  laid. 
Next,  a  hardened  steel  roller  was  mounted  to  the  milling 
machine  in  place  of  a  cutting  tool  so  that  it  acted  like  a 
rolling  pin.  The  process  was  similar  to  rolling  dough  with 
a  rolling  pin,  except  the  rolling  pin  spun  in  one  position 
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and  the  table  moved  back  and  forth.  Finally,  roll  fraction 
was  controlled  by  adjusting  the  table  height. 

The  method  of  actually  rolling  the  screen  was  developed 
by  trial-and-error .  The  results  are  as  follows: 

1.  The  screen  was  cut  into  strips  approximately 
32  mm  X  248  mm  (1.25  inch  x  9.75  inch). 

2 .  One  strip  was  used  as  the  prototype  for  its 
mesh  and  roll  fraction  to  adjust  the  table  height. 

3.  Each  strip  was  flatted  with  four  passes  of  the 
roller.  Before  each  pass  it  was  flipped  either  end-for-end 
or  side-to-side.  Thus,  each  corner  of  the  strip  was 
eventually  in  each  corner  of  the  parallel  plate. 

4.  Each  strip  was  checked  for  thickness  with  a 
micrometer  in  six  locations  to  ensure  the  thickness  was 
within  ±  0.03  mm  (  ±  0.0001  inches). 

5.  After  all  the  strips  of  like  mesh  and  percent 
roll  were  prepared,  they  were  stored  in  a  large,  marked 
envelope  until  further  use. 

A  discussion  of  the  effect  rolling  had  on  the  screen 
with  photographs  is  included  in  Chapter  IV. 

Screen  Cutting.  Cutting  the  screen  was  a 
relatively  simple  process,  compared  to  rolling  it.  The 
process  required  a  hardened  steel  punch,  a  hammer,  and  a 
smoothly  planed  hard  maple  board.  Three  strips  of  rolled 
screen  were  laid  on  top  of  each  other  then  taped  to  the 
board.  The  punch  was  placed  squarely  on  top  of  the  screen 
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and  all  three  strips  were  punched  simultaneously  with  a 
hammer.  Approximately  fifty  disks  were  allowed  to  back  into 
the  tube  portion  of  the  punch  and  then  the  stack  was  pushed 
out  of  the  tube  into  an  appropriately  marked  bottle.  The 
dimensions  of  the  punch  are  shown  in  Fig  3.4. 


Figure  3.4. 


Dimensions  of  Cylindrical  Punch  Used  to  Cut 
Screen  Disks 


The  quantity  and  quality  of  the  screen  disks  were 
closely  monitored  during  the  entire  process  of  regenerator 
assembly.  Any  unsatisfactory  disks  were  set  aside  and  were 
subtracted  from  the  total  count  of  disks  cut  (number  of 
holes  X  number  of  layers).  Unsatisfactory  disks  were  usually 
the  result  of  incomplete  cutting. 

Screen  Cleaning.  After  the  screens  had  been 
rolled  and  cut,  machine  oils  and  fingerprints  were  removed. 
All  of  the  screens  of  a  particular  mesh  and  percentage  roll 
were  placed  into  a  glass  jar  and  were  completely  covered 
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with  acetone.  Next,  the  entire  contents  was  swirled  for  at 
least  15  seconds  and  the  acetone  discarded.  Then,  a  second 
batch  of  clean  acetone  was  added  to  the  still  "wet"  screen 
disks  and  the  swirling  was  repeated.  After  the  second 
rinse,  the  acetone  was  discarded  and  the  screen  disks  were 
piled  on  a  heavy  duty  paper  shop  towel  to  air  dry.  Once 
dry,  latex  gloves  were  used  to  return  the  screens  to  their 
marked  bottle,  and  the  combined  mass  of  the  bottle  and 
screen  disks  was  recorded. 

Regenerator  Assembly .  Finally,  with  the  screen 
prepared,  the  entire  assembly  was  ready  to  be  completed. 

The  final  assembly  process  spanned  at  least  three  days 
because  of  the  cure  times  required  for  the  epoxy.  Epoxy  was 
selected  over  lead  solder  to  bond  the  entire  assembly 
together.  This  was  because  of  the  difficulty  in  applying 
lead  solder  to  stainless  steel.  In  addition,  silver  solder 
was  ruled  out,  except  for  the  thermocouple  assemblies, 
because  of  potential  heat  distortion  to  the  thin  stainless 
steel  tiibe  regenerator  housing. 

As  shown  in  Fig.  3.3  the  outlet  end  of  the  regenerator 
contained  two  oversized  screens  which  were  cut  using  a 
standard  5/8  inch  (15.875  mm)  punch.  The  two  screens  held 
the  entire  stack  of  screen  disks  in  the  tube.  The  two 
oversized  screens  were  epoxyed  to  the  flat  surface  of  the 
brass  end  cap.  The  epoxy  cured  overnight  with  a  glue  clamp 
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firmly  pressing  the  tube  against  the  two  oversized  screen 
disks  and  brass  end  cap. 


Figure  3.5.  Photo  of  Exploded  View  of  Regenerator 


The  screen  disks  were  packed  into  the  tube  after  the 
epoxy  cured  overnight.  The  screen  disks  were  packed  in 
small  groups  of  10  to  20  by  hand  using  latex  gloves.  First, 
a  group  was  picked  up  and  lined  up  to  make  a  straight,  even 
stack.  The  entire  group  was  started  into  the  tube  end. 

Next,  a  push  rod  [14.6  mm  (1.575  inch)  dieuneter]  was  used  to 
press  the  group  against  the  previously  packed  screen  disks. 
After  several  groups  were  positioned,  force  was  applied 
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manually  to  the  rod  to  ensure  the  discs  were  tightly  packed. 
This  entire  process  was  continued  until  no  more  disks  could 
be  placed  into  the  tube.  The  distance  from  the  top  screen 
to  the  top  of  the  tube  was  measured  and  recorded. 

The  unused  screen  disks  were  counted  and  returned  to 
their  designated  bottle.  Thus,  the  number  of  screen  disks 
packed  in  each  regenerator  was  calculated  from  the  count 
made  while  cutting  the  screen.  In  addition,  the  mass  of  the 
bottle  with  the  remaining  screen  disks  was  measured,  so  the 
mass  of  screens  in  the  regenerator  could  be  calculated.  The 
accuracy  was  better  than  0.04  percent. 

Finally,  the  regenerator  was  ready  for  the  final  step 
in  the  assembly  process.  The  inlet  end  cap  was  epoxyed  into 
place,  axial  pressure  was  applied  with  a  glue  clamp  and  the 
epoxy  was  allowed  to  cure  overnight. 

Pressure  Test.  The  final  step  in  the  regenerator 
assembly  was  to  perform  a  pressure  test.  This  ensured  the 
safety  of  the  regenerator,  under  the  testing  pressures  of  30 
atmospheres.  The  pressure  to  the  assembly  was  gradually 
increased  to  34  atmospheres  (500  psig) .  Then,  the  assembly 
was  checked  for  leaks  using  a  soap  suds  solution.  Every 
regenerator  assembly  passed  this  test  and  no  leaks  were 
found . 
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Flow  Meter 

The  flow  meter  used  to  measure  the  mass  flow  rate  of 
the  helium  gas  was  a  Brooks  Full  View  Rotameter  model 
9-1112-10,  serial  number  6412-75276.  The  meter  was 
calibrated  in  standard  cubic  feet  per  minute  (SCFM)  air. 

The  manufacture's  calibration  plate  stated  that  the 
SCFM^j^f  =  0.1345  X  %scale .  In  addition,  standard  atmosphere 
was  defined  as  70  ^F  and  0  psig.  Equation  2.22  provided  the 
conversion  to  kg/s  of  helium  from  the  rotameter  scale.  The 
minimum  scale  reading  was  10  percent  scale  which 
corresponded  to  approximately  0.28xl0“3  kg/s  and  the 
graduations  were  every  2.5  percent.  The  full  scale  reading 
of  100  percent  corresponded  to  approximately  3.2x10“"^  kg/s. 

Pressure  Transducers 

The  system  contained  three  different  pressure 
transducers  all  of  which  supplied  a  voltage  to  a  digital 
voltmeter  so  that  the  readings  could  be  recorded.  The  three 
transducers  measured  inlet  pressure,  Pj}  pressure  drop 
across  the  regenerator,  AP;  and  flow  meter  pressure, 

The  next  three  paragraphs  state  the  range,  transducer 
manufacturer,  model,  signal  conditioner,  and  digital 
voltmeter  used  for  each  of  the  three  pressure  measurements. 
In  addition,  the  calibration  technique  and  equipment  are 
described.  The  actual  calculations  made  with  these  pressure 


3-12 


measurements  are  included  in  Chapter  II,  Theory  and 
Experimental  Analysis. 

Inlet  Pressure.  P^.  The  Pj  transducer  was  required  to 
measure  the  maximum  inlet  pressure  of  30  atm,  or  roughly  450 
psia.  A  differential  pressure  was  used  instead  of  an 
absolute  pressure  transducer  because  of  cost  and  delivery 
time.  In  addition,  a  digital  barometer  was  readily 
available  so  that  gage  pressure  readings  could  be  easily 
converted  to  absolute  pressure.  A  Validyne  DP15-56,  0-500 
psid  pressure  transducer,  serial  number  84654,  was  selected 
for  its  accuracy  and  reasonable  cost.  It  was  connected  to  a 
Validyne  Dual  Modulator,  CD  280,  serial  niunber  109040.  The 
voltage  from  the  CD  280  was  measured  with  a  Hewlett  Packard 
HP  3466A  Digital  Multimeter,  serial  number  1716A-18320. 

Pressure  Drop,  AP.  A  CEC  model  4-351-002,  0  to  50  psid 
pressure  transducer,  serial  number  2357,  was  selected  to 
measure  the  pressure  drop  through  the  regenerator.  Its 
range  was  satisfactory  and  it  was  "on-hand".  In  the  process 
of  testing,  two  transducers  failed  so  a  total  of  three 
identical  transducers  were  used  (the  calibration  for  each 
transducer  was  closely  tracked) .  The  second  and  third 
serial  numbers  were  2683  and  2769,  respectively.  In  all 
three  cases,  the  power  was  supplied  with  an  Endevco  Power 
Supply,  model  4225,  serial  number  AG  23;  and  the  signal  was 
amplified  with  an  Endevco  Signal  Conditioner,  model  4423, 
serial  number  BA  65.  The  voltage  was  displayed  on  a  Hewlett 
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Packard  HP  3466A  Digital  Multimeter,  serial  number  1716A- 
18308. 

Flow  meter  Pressure.  The  pressure  of  the  helium 

entering  the  flow  meter  was  expected  to  be  slightly  above 
atmospheric  pressure  because  of  back  pressure  from  the 
exhaust  tube  line.  Once  again,  barometric  pressure  was 
readily  available  so  a  differential  pressure  transducer 
configured  for  gage  pressure  was  used.  A  Validyne 
Differential  Pressure  Transducer,  0-20  psid,  model  DP-10-42, 
serial  number  74960  was  selected.  Its  range  was  accept6ible 
and  it  was  compatible  with  the  Validyne  Dual  Modulator, 

CD  280,  serial  niimber  109040  (same  modulator  as  Pj,  but 
different  channel).  The  voltage  was  displayed  on  a  TSI 
Voltmeter,  Model  1076. 

Pressure  Calibration.  All  five  pressure  transducers 
were  calibrated  using  dead  weight  testers.  All  the 
associated  equipment  listed  above  for  each  transducer  was 
used  during  the  calibration.  The  voltage  vs.  pressure  was 
measured  for  a  total  of  more  than  11  points  in  the  range  of 
interest.  The  readings  were  taken  as  the  pressure  increased 
then  decreased.  An  Ametek,  Type  K  Pneumatic  Pressure 
Tester,  model  HK-500,  serial  number  79672,  used  for  the  Pj 
and  AP  transducers.  An  Ametek,  Type  K,  Pneumatic  Dead 
Weight  Pressure  Tester  model  PK2-254WC,  serial  number  79055 
was  used  for  the  P^^^  transducer. 
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Barometric  Pressure.  Barometric  pressure  was  recorded 
from  a  CEC  Transamerica  Delaval  Digital  Barometer  model 
2500-0103,  serial  number  258.  The  barometer  display  was  in 
millimeters  of  mercuiry. 

Temperature 

Temperature  measurement  at  the  regenerator  inlet  and 
outlet  were  a  crucial  part  of  this  study.  The  next  few 
paragraphs  explain  the  selection  of  type  T  thermocouples  and 
the  associated  equipment  which  processed  and  stored  the 
data.  In  addition,  the  software  which  operated  the 
equipment  is  explained. 

Selection  of  Type  T  Thermocouples.  Temperature 
measurements  needed  to  be  fast  and  accurate  to  measure  the 
transient  temperature  response  which  resulted  from  the  step 
temperature  change  (7^-f,).  Thermocouples  provided  the 
millisecond  response  required.  Type  T  (copper-constantan) 
thermocouples  were  acceptable  for  the  temperature  range  used 
in  this  experiment.  The  thermocouples  used  were  Omega 
Subminiature  Quick  Disconnect  Style  Sheathed  Thermocouple 
Probes  with  exposed  junctions,  model  TMQSS-032(Z)-6.  The 
sheath  diameter  was  0.032  inches  (0.813  mm)  and  the  wire 
diameter  was  0.006  inches  (0.145  mm).  The  sheath  and  wire 
diameters  were  adequate  enough  to  allow  them  to  resist  the 
drag  affects  from  the  flow.  However,  the  junction  and  wire 
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diameters  were  small  enough  to  provide  the  required  quick 
response . 

Figure  3.6  illustrates  the  thermocouple  assembly.  The 
inlet  and  outlet  assemblies  were  identical.  The  inlet  and 
outlet  thermocouples  were  always  used  in  the  same  position. 

A  third  temperature,  Tj,  was  also  measured.  A  matching 
subminiature  thermocouple  was  epoxyed  into  a  small  hole  in  a 
Swagelok  3/8  inch  union  fitting  with  "Double  Bubble"  brand 
epoxy.  The  fitting  was  approximately  120  mm  from  the  four¬ 
way  valve  in  the  secondary  flow. 


1/4  in.  (6.35  mm) 
Stainless  Steel 
Tube 


Silver 

Solde 


Quick  Disconnect 
Thermocouple 
Type  T 

—  Thermocouple 
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Thermocouple 
Junction 


6.00  in.  (15.3  mm) 

Modified 
No.  10  Bolt 

0-Ring 

0.46  in.  (12  mm) 


Figure  3.6.  Thermocouple  Assembly 


Signal  Processing.  Three  pieces  of  equipment  were 
required  to  process  the  thermocouple  output  into 
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temperature.  First,  a  Kelthley  MetraByte  Corporation  EXP-16 
16-channel  analog  multiplexer  (serial  niimber  233266) 
amplified  and  multiplexed  the  analog  signal  for  the  analog 
to  digital  (A/D)  converter.  The  A/D  converter  was  a 
Keithley  MetraByte  Corporation  DAS-8  12  bit  A/D  converter 
(serial  number  9537).  The  DAS-8  was  connected  to  a  Zenith 
Z-248  International  Business  Machines  (IBM)  compatible 
computer . 

The  manufacturer  recommended  gain  setting,  for  type  T 
thermocouples,  was  200.  This  included  the  entire  normal 
operating  range  of  type  T  thermocouples.  The  gain  resulted 
in  a  temperature  resolution  of  0.3022  K.  However,  this 
study  was  for  a  much  smaller  range  of  temperatures.  Thus, 
the  maximum  gain  of  1000  was  used  to  improve  the  resolution 
to  0.0608  K. 

Temperature  Acquisition  Software.  The  software  used  to 
collect  the  temperature  data  was  written  in  6W-BASIC  3.20. 
The  manufacturer  provided  sample  software  for  the  user  to 
modify  for  the  user's  needs.  The  program  titled  "EXP-T.BAS" 
was  the  core  program  used  to  build  the  software  used  for 
temperature  acquisition.  This  program  performed  the 
functions  necessary  to  retrieve  temperature  data  from  up  to 
16  temperature  ports  in  the  EXP-16.  The  program  also 
displayed  these  temperatures  on  the  computer  monitor  in  real 
time.  The  author  streamlined  and  modified  the  sample 
software  so  that  it  met  the  specific  needs  of  this  study. 


3-17 


The  final  program  (listed  in  Appendix  E)  had  three  modes  of 
operation.  The  three  modes  were  file  initiation,  real  time 
display  and  data  acc[uisition. 

The  first  mode,  file  initiation,  allowed  the  user  to 
enter  the  filename  of  the  test  which  was  to  be  run.  The 
program  initiated  the  matching  filenames  on  the  a,  b,  and  c 
drives  of  the  computer.  The  program  then  returned  the 
computer  to  real  time  display  mode. 

The  real  time  display  mode  displayed  the  three 
thermocouple  temperatures,  the  cold  junction  temperature, 
the  filename,  and  status  of  the  current  file  ("old"  or 
"new").  This  was  the  default  mode  whenever  one  of  the  other 
two  modes  were  not  in  progress.  The  other  two  modes  were 
entered  from  the  real  time  display  mode.  This  mode  allowed 
the  user  to  know  when  conditions  were  favorable  to  perform  a 
test.  It  continuously  updated  and  displayed  the 
temperatures  until  the  file  initiation  or  data  acquisition 
mode  was  selected. 

The  data  acquisition  mode  was  streamlined  to  maximize 
the  rate  of  data  acquisition.  The  display  of  temperatures 
was  no  longer  updated  until  an  entire  test  was  completed. 

The  following  sequence  of  events  occurred.  First,  the  cold 
junction  temperature  was  recorded  (experience  showed  the 
cold  junction  temperature  to  be  steady  during  the  18  seconds 
required  to  collect  data) .  Second,  the  program  entered  a 
loop  which  stored  the  clock  time  and  three  thermocouple 
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voltages  for  each  of  the  500  readings.  Third,  the  program 
converted  the  clock  time  to  elapsed  seconds  and  converted 
thermocouple  voltages  to  temperatures.  The  program  saved 
these  calculated  values  to  the  data  file  on  each  disk  drive. 
Finally,  the  program  returned  to  real  time  display  mode. 

Tubing 

The  factors  used  to  select  the  type  and  size  of  tubing 
were  operating  pressure,  flow  velocity,  availability,  cost 
and  workeJsility.  The  tubing  had  to  be  pressure  rated  to  at 
least  30  atm  (3  MPa)  and  have  a  large  enough  inside  diameter 
to  keep  flow  velocities  low  enough  to  eliminate  Mach  number 
compressibility  effects  (actual  maximum  velocities  were 
approximately  100  m/s  or  Mach  number  equal  to  0.1).  The 
considerations  of  availability,  cost  and  workability  all 
concluded  in  using  copper  tubing.  Thus,  the  final  selection 
was  3/8  inch  (9.525  mm)  outside  diameter  with  0.032  inch 
(0.81  mm)  wall  thickness  soft  copper  tubing. 

Valves 

The  four  types  of  valves  used  in  the  test  apparatus 
were  a  four-way  valve,  3  two-way  valves,  2  metering  valves, 
and  a  pressure  regulator.  All  the  valves  required  a 
pressure  rating  of  at  least  30  atm  (3  MPa)  and  had  to  be 
sized  to  allow  sufficient  flow  rates.  The  sizing  criteria 
was  flow  coefficient,  Q.  The  methods  of  calculating  Q 
vary  slightly  between  manufacturers.  As  a  results,  care  was 
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required  when  selecting  each  valve.  The  next  few  paragraphs 
describe  the  basic  conditions  for  which  each  valve  was 
selected.  The  temperature  used  to  find  each  was  70°  F 
(294  K)  but  the  pressure  was  the  lowest  expected  at  each 
valve's  location. 

Four-Wav  Valve.  The  four-way  valve  had  to  quickly 
switch  the  flow  to  provide  the  step  temperature  change 
required  for  the  experiment.  Each  channel  was  sized  to 
allow  approximately  3  x  10”^  kg/s  flow  at  an  outlet  pressure 
of  10  atm.  The  valve  selected  was  a  Whitey  4-Way  Ball 
Valve,  model  B-45yF8. 

Two-Way  Valves.  The  3  two  way  valves  allowed  quick 
shut-off  and  isolation  of  the  system  to  preserve  helium. 

The  two-way  valves  were  sized  to  allow  approximately 
3  x  10“3  kg/s  flow  at  10  atm.  The  valve  selected  was  a 
Whitey  2-Way  Ball  Valve,  model  B-44S6. 

Metering  Valves.  There  were  two  different  metering 
valves.  Both  were  required  to  withstand  30  atm  differential 
pressure.  The  was  calculated  for  3  x  10”3  kg/s  flow  at 
6.5  atm  (to  allow  for  pressure  drop  experienced  upstream  in 
the  regenerator) .  The  toughest  requirement  was  to  have  a 
smooth,  linear  response  of  versus  n\unber  of  turns  so  that 
the  flow  rate  could  be  precisely  controlled.  The  primary 
flow  was  controlled  with  a  Nupro  Metering  Valve,  model  SS- 
4BMRG.  The  secondary  flow  metering  valve  had  to  meet  many 
of  the  s£une  requirements  as  the  primary  valve  but  was  less 
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expensive.  The  secondary  flow  was  controlled  with  a  Whitey 
Screwed  Needle  Valve,  model  B-4R56. 

Pressure  Regulator.  The  pressure  regulator  needed  to 
reduce  the  tank  pressure  of  2300  psi  (16  MPa)  to  the  10  to 
30  atm  (150  to  450  psi)  range  required  for  the  tests.  In 
addition,  the  valve  required  a  large  enough  orifice  to  allow 
the  required  total  flow  rate  of  approximately  3  x  10~^  kg/s. 
A  Grove  15-L,  3000  psi  inlet  with  0-750  psi  outlet,  serial 
number  L-41575  was  used. 

Ice  Bath  Design 

The  cold  temperature  flow  was  created  by  chilling  room 
temperature  helium  with  an  ice  bath.  The  heat  exchanger 
surface  was  made  from  approximately  3.2  m  of  tubing  rolled 
into  a  0.17  m  diameter  coil.  The  tubing  was  the  same  type 
of  copper  tubing  used  in  the  rest  of  the  apparatus.  The 
entire  device  wan  housed  in  a  Playmate  18  can  cooler.  The 
cooler  was  filled  with  ice  provided  free  of  charge  by  the 
local  snack  bar. 

The  ice  bath  design  was  based  on  several  simplifying 
assumptions  to  estimate  the  required  txibing  length.  The 
first  assumption  was  a  constant  wall  temperature  of  0 
(273  K) .  The  second  assumption  was  a  straight  tube  with 
laminar  flow  (worst  case).  Finally,  an  entrance  temperature 
of  300  K  and  an  outlet  temperature  of  274  K  was  assumed. 
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Unfortunately,  the  temperatures  entering  the  regenerator 
were  between  285  K  and  290  K. 

Helium 

The  helium  was  "high  purity  helium"  manufactured  by  the 
United  States  Bureau  of  Mines,  Helium  Field  Operations.  The 
heliim  was  at  least  99.997  percent  pure.  The  helium  met  the 
requirements  of  MIL-P-27407A  (Type  I,  Grade  A). 

Insulation 

Insulation  was  added  to  insulate  the  line  between  the 
ice  bath  and  the  four-way  valve.  In  addition,  a  removable 
section  of  insulation  was  placed  over  the  regenerator  which 
extended  to  the  four-way  valve.  The  insulation  was  standard 
black  foam  used  for  pipe  insulation.  The  foam  tube  was  3/8 
in.  (9.5  mm)  thick.  The  insulation  was  tightly  wrapped 
around  the  tubing  and  secured  with  wire  bundle  ties  to  form 
a  seal. 

Fittings 

All  fittings  were  required  to  be  safe  for  the  operating 
pressure  of  30  atm  and  to  fit  the  3/8  inch  (9.525  mm) 
outside  diameter  copper  tubing.  In  addition,  they  needed  to 
be  able  to  survive  numerous  assemblies  and  disassembles  when 
helium  tanks  and  regenerators  were  changed.  Finally,  they 
were  required  to  accommodate  the  various  valves  and  other 
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connections  in  the  system.  The  final  selection  was  Swagelok 
brand  brass  fittings. 

Leak  Checks 

All  fittings  and  connections  were  checked  for  leaks 
upon  initial  assembly  and  whenever  they  were  opened  or 
altered  in  any  way.  All  leak  checks  were  perfonned  by 
brushing  a  small  amount  of  soap  suds  solution  over  the  area 
in  question.  Any  bubbles  or  foaming  was  an  indication  of  a 
leak  and  corrective  action  was  taken. 
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IV.  Experimental  Results 


The  previous  chapters  explained  the  background,  theory 
and  experimental  design.  These  chapters  lead  to  the 
experimental  results.  This  chapter  is  divided  into  three 
main  sections.  The  first  topic  is  the  physical  effect  that 
rolling  had  on  the  screen.  The  second  topic  is  the  effect 
rolled  screen  had  on  friction.  The  third  topic  is  the 
effect  rolled  screen  had  on  heat  transfer.  Finally,  the 
combined  effects  of  friction  and  heat  transfer  are 
discussed. 

In  order  to  simplify  the  discussion,  the  fraction  roll, 
R,  is  discussed  in  terms  of  the  approximate  values  of  0.00, 
0.15,  0.30,  and  0.50.  In  actuality,  the  real  values  were 
slightly  different  for  250  and  325  mesh  screens.  This 
simplification  eases  the  comparison  between  250  and  325  mesh 
screen.  However,  all  graphs  refer  to  the  actual  measured 
value  of  R. 

Physical  Effect  of  Rolling 

The  physical  effect  rolling  had  on  screen  is 
illustrated  by  scanning  electron  microscope  photographs  in 
Figs.  4.1  and  4.2  (250  and  325  mesh  respectively).  The 
first  difference  between  the  two  figures  is  that  the  250 
mesh  was  plain  weave  and  the  325  mesh  was  twill  weave.  This 
was  because  of  the  unavailability  of  plain  weave  325  mesh 
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screen.  However,  both  sets  of  photographs  display  the  same 
effect  due  to  rolling. 

For  both  the  250  and  325  mesh  unrolled  (0  fraction 
roll)  screens  the  cross  section  of  the  individual  wires  was 
nearly  circular  in  all  locations.  The  effect  of  rolling  was 
barely  evident  with  the  0.15  fraction  roll  screen.  The 
distortion  was  concentrated  on  the  top  of  each  weave 
resulting  in  a  slight  flat  spot.  However,  neither  the 
contact  region  between  crossing  wires  or  the  square  gap 
between  each  weave  showed  any  noticeable  change.  The  0.30 
fraction  roll  screen  demonstrated  the  same  effects  as  the 
0.15  fraction  roll  except  that  the  flattened  region  at  the 
top  of  each  weave  was  larger. 

Tne  0.50  fraction  roll  screen  showed  the  first  signs  of 
large  physical  distortions  due  to  rolling.  First,  the 
flattened  region  at  the  top  of  the  weave  was  significantly 
larger.  Second,  there  was  some  distortion  in  the  contact 
region  between  crossing  wires.  Finally,  the  square  gap 
between  each  weave  was  noticeably  reduced  and  no  longer  had 
square  corners.  Thus,  the  effect  of  rolling  the  screen  was 
small  for  the  0.15  and  0.30  fraction  roll  screen,  but  showed 
significant  effect  for  0.50  fraction  roll. 

Friction  Results 

This  section  covers  the  results  of  the  friction 
measurements.  First,  the  graphical  result  of  friction 


factor,  /,  versus  Reynolds  nvunber  for  internal  flow.  Re,  is 
discussed.  Second,  the  repeatability  results  are  presented. 
Third,  a  method  to  predict  friction  factor  of  rolled  screen 
is  presented.  Fourth,  the  variation  of  drag  coefficient, 

C[j,  with  Reynolds  number  for  external  flow,  Re^,  is 
presented  to  include  correlation  with  previous  results. 

Then,  tests  from  a  hollow  regenerator  are  presented. 

Finally,  the  overall  effect  of  rolled  screen  is  discussed. 

Friction  Factor,  f.  The  standard  method  to  present 
friction  results  in  Kays  and  London  (1984)  was  the  friction 
factor,  /.  The  investigator  used  the  method  in  Chapter  II 
to  convert  the  experimental  test  data  to  /  vs.  Re.  The  next 
two  paragraphs  and  figures  demonstrate  sample  results  from 
tests  with  "unchoked"  and  "choked"  flow.  Then  the  combined 
results  showing  trends  from  rolling  are  presented. 

Figure  4.3  illustrates  the  friction  results  for  the 
tests  of  250  mesh  screen  with  R  =  0.15.  This  figure  is  an 
example  of  tests  performed  when  none  of  the  tests 
experienced  "choked"  flow  (choked  flow  is  explained  in  the 
next  paragraph).  As  Fig.  4.3  illustrates,  /  is  not 
influenced  by  variations  in  inlet  pressure.  The  slight 
divergence  near  Re  =  100  was  the  result  of  random  noise 
since  the  pattern  is  not  consistent  with  any  of  the  other 
series  of  tests. 

Figure  4.4  illustrates  the  friction  results  for  the 
tests  of  325  mesh  and  R  »  0.50.  This  figure  ic  an  example 
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Figure  4.3.  Friction  Factor,  /,  vs.  Reynolds  NusUoer,  Re, 
for  250  Mesh  Screen  with  0.143  Fraction  Roll.  Uncertainty 
ranged  from  12  percent  at  the  lowest  Re  and  decreased  to  1.1 
percent  at  the  highest  Re  for  each  curve. 
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Figure  4.4.  Friction  Factor,  /,  vs.  Reynolds  Nus^r,  Re, 
for  325  Mesh  Screen  with  0.485  Fraction  Roll.  Uncertainty 
ranged  from  12  percent  at  the  lowest  Re  and  decreased  to  1.1 
percent  at  the  highest  Re  for  each  curve. 


4-6 


of  the  worst  case  for  choked  flow.  Choked  flow  occurred 
«rhen  the  pressure  drop  across  the  regenerator  was  such  that 
the  flow  rate  could  no  longer  be  increased  by  opening  the 
flow  control  valve.  The  pressures  exiting  the  regenerator 
in  all  of  the  choked  flow  cases  was  greater  than  3  atm  and 
the  systesi  vented  to  1  ata.  This  exit  pressure  of  3  atm 
suggests  that  the  flow  was  not  choked  by  the  regenerator 
Itself.  Thus,  the  flow  was  choked  by  the  smallest  orifice 
in  the  downstream  system  which  was  the  flow  control  valve 
body.  The  teats  illustrated  in  Fig.  4.4  experienced  choked 
flow  at  F;  «  10  and  20  atm.  The  regenerators  which 
experienced  choked  flow  were  250  mesh  R  •  0.50;  and  325  sMsh 
R  -  0.15,  0.30,  0.50.  These  regenerators  experienced  choked 
flow  trhen  the  inlet  pressure  was  10  ata.  The  325  mesh,  R  > 
0.50  regenerator  also  choked  at  20  ata.  In  all  choked 
cases,  /  was  loiser  than  the  unchoked  friction  factors. 

The  trend  caused  by  rolling  the  screen  is  shown  in 
figs.  4.5  and  4.6.  In  both  figures  the  curves  shifted  to 
the  left  and  down  as  the  roll  fraction  increased.  The  two 
basic  trends  were  the  shift  in  both  Re  and  /. 

As  shown  in  Chapter  II,  Rc  was  a  function  of  both  mass 
flow  rate,  ff,  and  hydraulic  radius,  The  mass  flow  rates 

resMiined  roughly  the  sasM  for  each  grouping  of  data.  For 
exaa4>le,  all  regenerators  used  •  0.28  x  I0~^  kg/s  and 

•2.9  X  10“3  kg/s.  So,  the  shift  to  the  loft  on  Pigs. 
4.5  and  4.6  was  due  to  the  decrease  in  caused  by  rolling. 
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Figure  4.5.  Friction  Factor,  f,  vs.  Reynolds  Niunber,  Re, 
for  250  Mesh  Screen.  Uncertainty  ranged  from  12  percent  at 
the  lo%»e8t  Re  and  decreased  to  1.1  percent  at  the  highest  Re 

for  each  curve. 


Figure  4.6.  Friction  Factor,  /,  vs.  Reynolds  Number,  Re, 
for  325  Mesh  Screen.  Uncertainty  ranged  from  12  percent  at 
the  lowest  Re  and  decreased  to  1.1  percent  at  the  highest  Re 

for  each  curve. 


The  decrease  in  f  between  regenerators  can  be 
misleading.  The  pressure  drop  across  two  different 
regenerators  with  the  same  outside  dimensions  is  not  only  a 
function  of  f,  but  also  the  porosity,  p,  and  the  total  area, 
A.  A  method  to  estimate  /  and  Re  is  explained  in  a  later 
section . 

Repeatability.  The  repeatability  is  illustrated  in 
Fig.  4.7.  The  repeatability  tests  were  performed  using  the 
250  mesh  R  =  0.30  regenerator  at  -  20  atm.  The 
regenerator  was  connected  in  the  same  manner  as  all  other 
tests.  The  original  data  were  collected  on  the  same  day  but 
the  repeat  data  were  collected  on  two  different  days.  In 
addition,  for  each  different  day's  data  the  250  mesh 
R  -  0.30  regenerator  was  freshly  reinstalled.  As  shown  in 
Fig.  4.7,  the  repeatability  was  nearly  perfect. 

Correction  Factor  for  Rolling.  The  nearly  parallel 
lines  in  Figs.  4.5  and  4.6  implied  that  correction  of 
unrolled  screen  friction  factor  for  rolled  screen  was 
possible.  The  investigator  defined  the  subscript  "0"  as  the 
property  of  unrolled  screen  and  no  subscript  as  the  property 
for  rolled  screen.  This  resulted  in  the  following 
definitions: 

/  *  friction  factor 

»  friction  factor  for  unrolled  screen 

^  *  Reynolds  number  for  internal  flow  at  the  same 
mass  flow  rate  as  Re^ 

Re^x  Reynolds  number  for  Internal  flow  for  unrolled 
screen 
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p  *  porosity 

Pg  -  porosity  of  unrolled  screen 


Figure  4.7.  Repeatability  Test  With  250  mesh  0.286  Fraction 
Roll.  Friction  Factor,  /,  vs.  Reynolds  Number,  Re. 
Uncertainty  ranged  from  12  percent  at  the  lowest  Re  and 
decreased  to  1.1  percent  at  the  highest  Re  for  each  curve . 


The  investigator  developed  a  method  to  normalize  the 
results  in  Figs.  4.5  and  4.6.  The  resulting  normalized 
friction  factor  was  fp/f„p„.  Figure  4.8  shows  that  fp/foPo  vs. 
fraction  roll,  R,  converged  into  a  nearly  linear 
relationship.  Because  the  lines  of  constant  Re^  overlaped, 
the  lines  were  removed  and  only  the  data  points  were 
plotted. 

The  six  circled  points  in  Fig.  4.8  ware  all  at 
P,  "  10  atm  and  at  the  highest  flow  rate  tested.  Four  of 
the  six  circled  points  were  at  the  choked  flow  conditions 
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mentioned  previously  and  the  other  two  were  measured  at  the 
maximum  flow  rate  of  2.9  x  10~3  kg/s. 


Figure  4.8.  Normalised  Friction  Factor  vs.  Fraction  Roll,  R 


Figure  4.8  depicts  two  different  linear  regressions. 

The  first  line  was  with  the  six  circled  points  and  the 
second  was  without  the  six  circled  points.  The  regression 
was  perfonsed  so  that  the  Intercept  of  the  vertical  axis  was 
forced  to  equal  one.  The  resulting  equations  were 

-^  =  l-0  845R,  all  points  (4.1) 

foP. 

=  ]-0904R ,  without  circled  points  (4.2) 

foPo 
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Rat-^matlon  Technicme  for  f.  For  250  and  325  mesh,  the 
investigator  developed  the  following  process  to  calculate  /: 
First,  calculate  Re^  from  properties  of  unrolled  screen 
using  Eq.  2.23.  Second,  calculate  /„  by  multiplying  C^, 
(improved  curve  fit  ,  Eq.  4.5,  or  some  other  source)  by //C^, 
(Eq.  2.5).  When  calculating  ^  use  Re^  and  p^.  Third,  either 
measure  p  if  the  regenerator  is  already  constructed  or 
approximate  it  using  p  •  p^{\-^R)-R  [derived  from  definitions 
of  p  (Eq.  2.15)  and  R  (Eq.  2.18)  assuming  perfect 
regenerator  packing].  Fourth,  apply  Eq.  4.1  or  4.2  to 
calculate  /.  Finally,  if  Re  is  desired  it  can  be 
approximated  as  Re  •  Rc„{l-R)  assuming  that  the  screen  can  be 
perfectly  and  evenly  packed. 

This  technique  is  useful  for  calculating  /  for  the 
parameters  of  this  research  (250/325  mesh,  10  <  Re^  <  100 
and  R  <  0.50  for  helium  gas).  However,  the  investigator 
expects  this  technique  could  be  extrapolated  to  other  mesh 
sises,  gases,  and  Reynolds  nusO^ers. 

Comparison  to  Previous  Research.  As  discussed  in 
Chapter  II,  Tong  and  London  (1957)  provided  the  following 
empirical  curve  fit  for  screen  matrixes: 


where 


Jog  10  Cq 


•  33  -0  3J 


0.54 

P 


(2.4) 


Cq  *  drag  coefficient  per  screen 

Re^  -  Reynolds  number  for  external  flow  (Eq.  2.24) 
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Legend 


-f-  Data  250  mesh 


Figure  4.9.  Ccunparlson  to  the  Curve  Fit  of  Tong  and  London 
(1957).  Drag  coefficient,  Q,,  vs.  Reynolds  Number.  Rc^. 
Uncertainty  ranged  from  12  percent  at  the  lo%rast  Re^  and 
decreased  to  1.1  percent  at  the  highest  Re^  for  each  curve. 


The  Tong  and  London  (1957)  curve  fit  is  plotted  on  Fig. 
4.9  with  the  data  for  unrolled  screen.  The  agreement  was 
good.  The  Tong  and  London  (1957)  data  was  taken  with  air 
using  a  significantly  larger  matrix  of  crossed  rods  and 
corrected  for  woven  screen. 

A  correction  to  the  curve  fit  of  Tong  and  London  (1957) 
was  applied  by  using  an  equation  of  the  form 

^10^0  =  — f-RCrf*'  +  — ^  .  (4*3) 

P  P 

The  constants  A^,  and  were  to  be  found.  A  standard 
least  squares  technique  (Strang,  1988tl54)  could  not  be 
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applied  to  solve  for  all  three  constants  simultaneously. 

This  was  because  the  right  hand  side  of  Eq.  4.3  could  not  be 
separated  into  linearly  independent  terms.  The  investigator 
attempted  several  iteration  techniques  but  the  best  solution 
was  found  by  assuming  Tong  and  London's  value  of  =  -0.54. 
This  resulted  in  a  linear  regression  of 

lo8io(p^  logio Cd-pCo)  =  logjo  iog,o  Re^  .  (4.4) 
Thus,  the  linear  regression  of  Eq.  4.4  using  unrolled  screen 
data  provided  the  new  curve  fit 


logic 


P 


0.54 

P 


(4.5) 


As  is  shovm  in  the  next  section,  this  relation  holds  for 
R  0.30  using  p  =  p^.  The  new  curve  fit  is  also  plotted  on 
Fig.  4.9. 

Drag  Coefficient  Results.  The  drag  coefficient,  C^, 
was  calculated  by  converting  the  friction  factor,  /,  to 
using  Eq.  2.5.  In  addition,  the  Reynolds  number  for 
external  flow,  Re^,  was  calculated  using  Eq.  2.24.  Re^ 
remained  constant  for  a  particular  flow  rate,  Wf  and  mesh 
size.  Figs.  4.10  and  4.11  are  the  resulting  graphs. 

The  250  mesh  results  in  Fig.  4.10  are  quite 
interesting.  For  R  up  to  0.30  the  drag  coefficient  was 
unchanged.  In  addition,  for  325  mesh  screen  (Fig.  4.11)  Cq 
was  unchanged  for  R  *  0.15  and  it  increased  slightly  for 
R  •  0.30.  This  means  that  for  the  two  mesh  sizes  considered 
the  screen  can  be  rolled  up  to  R  >  0.30  without  any 
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Fraction  Roll 


Figure  4.10.  Drag  coefficient,  Q,,  vs.  Reynolds  Number, 
Re^,  for  250  mesh  screen.  Uncertainty  ranged  from  12 
percent  at  the  lowest  Re^  and  decreased  to  1.1  percent  at 
the  highest  Re^  for  each  curve. 


Figure  4.11.  Drag  coefficient,  C^,,  vs.  Reynolds  Number, 
Re^,  for  325  mesh  screen.  Uncertainty  ranv_  jd  from  12 
percent  at  the  loi«est  Re^  and  decreased  to  1.1  percent  at 
the  highest  Re^  for  each  curve. 


noticeable  change  in  C^, .  Q,  is  discussed  further  in  the 
friction  discussion. 

Rmpty  Tube  Tests.  Tests  were  performed  on  an  empty 
tube  which  was  constructed  in  the  same  manner  as  the  ether 
regenerators.  The  only  difference  was  that  screen  disks 
were  not  added.  The  four  tests  were  performed  at  Py  =  10 
and  30  atm  and  W  =  0.28  x  10"^  and  2.9  x  10“^  kg/s.  Thus, 
the  empty  tube  tests  were  performed  at  the  minimum  and 
maximum  extremes  of  the  test  conditions.  The  pressure  drop 
was  compared  to  the  unrolled  250  mesh  results  at  matching  Pj 
and  W.  The  AP  was  between  0.2  and  0.4  percent  of  the  250 
mesh  results.  Thus,  the  design  of  the  regenerator  case  had 
a  minimal  effect  on  the  test  results. 

Friction  Discussion.  As  was  discussed  earlier,  the 
friction  factor  for  250  and  325  mesh  screen  can  now  be 
predicted  for  10  <  Re„  <  100  (Eg.  4.5  and  Eg.  2.18).  The 
new  Reynolds  n\imber  after  rolling  can  also  be  predicted. 

Re  *  Re0(l-/{).  However,  these  values  do  not  provide  a  good 
intuitive  relationship  for  predicting  the  pressure  drop  for 
a  regenerator  with  rolled  screen. 

Rcj  and  Cp  seem  to  be  superior  to  Re  and  /  for 
correlation.  The  advantage  to  Re^  is  that  it  is 
uninfluenced  by  porosity  where  Re  is  influenced.  Because 
is  constant  for  R  ^  0.30,  the  pressure  drop  across  a  given 
number  of  screens  will  be  constant  whether  the  screens  have 
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been  rolled  or  not.  Thus,  is  more  useful  than  /  which 
does  not  remain  constant. 

Heat  Transfer 

This  section  presents  the  heat  transfer  results  for  the 
research.  First,  the  graphical  results  of  Colburn  factor, 
jfj,  versus  Reynolds  number  for  internal  flow.  Re,  are 
discussed.  Then,  the  repeatability  results  are  presented. 
Next,  an  approximate  method  to  predict  the  Colburn  factor 
for  screen  is  presented  for  varied  amounts  of  rolling.  The 
results  for  tests  with  a  hollow  tube  are  also  discussed.  In 
addition,  the  actual  convection  heat  transfer  coefficient, 
h,  is  plotted.  Finally,  the  overall  results  are  discussed. 

Colburn  Factor.  /^.  Figure  4.12  is  an  example  of  the 
teat  results  for  a  regenerator.  The  test  results  were  for 
250  mesh  with  0.15  fraction  roll.  This  figure  shows  that 
the  heat  transfer  results  were  not  as  tightly  grouped  as  the 
friction  results.  In  addition,  the  curve  shows  was 
independent  of  inlet  pressure  because  the  results  were 
randomly  scattered  around  the  basic  trend.  Since  there  was 
no  noticeable  trend  for  tests  with  choked  flow,  there  is  no 
figure  showing  this  data. 

Figures  4.13  and  4.14  illustrate  the  overall  trend  of 
rolling  on  vs.  Re.  In  order  to  simplify  the  appearance  of 
the  graphs  the  data  points  %#ere  not  connected  as  the 
friction  graphs  were.  Instead,  the  solid  lines  represent  a 
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curve  fit  of  the  results  for  all  the  data  for  a  particular 
regenerator.  These  curve  fits  were  made  by  performing  a 
quadratic  curve  fit  to  the  logarithms  of  and  Re. 
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Figure  4.12.  Sample  Test  Results  Using  250  Mesh  0.143 
Fraction  Roll.  Colburn  Factor,  vs.  Reynolds  Mtimber,  Re, 
for  325  mesh  screen.  Uncertainty  ranged  from  15  percent  at 
the  lowest  Re  and  increased  to  50  percent  at  the  highest  Rc 

for  each  curve . 


Like  the  friction  results,  the  curves  shifted  down  and 
to  the  left  with  incrcjsinq  fraction  roil,  /V  As  explained 
in  the  friction  section,  the  shift  to  the  left  on  the  graph 
was  due  to  the  decreasing  hydraulic  radius,  redut  ;  r.q  Ri 

A  i»thod  to  predict  is  presented  ir  a  .ater  see  t  .  o  ' 
the  results  are  further  discussed  w.tt  'he  -a  i  i  r  ^  .  t  * 
convection  heat  transfer  ’♦*ffi 
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fitly  Th«  r«p*AtAbility  t«ata  for  h«at 
trAiiafar  taara  run  aiaultan«oualy  with  th«  friction 
r«pMt«bility  t«ata.  The  twata  war*  parfonaad  on  2^0  wtah 
acraan  with  0.30  fraction  roll  at  20  ata  inlat  praaaura 
Tha  original  taata  wara  conductad  on  tha  aaaw  day  but  tha 
rapaat  taata  wara  parforaad  on  two  ulffarant  daya  Tha 
raaulta  of  thaaa  taata  ara  illuatratad  in  Pig  4  15  Tha 
rapaatabi 1 i ty  waa  not  aa  outatandinq  aa  tha  frirtion  taata, 
twt  tha  raaulta  ahowad  raaaunabia  rapaatabi i i t y 
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haliuB  gas.  Finally,  the  range  of  Re  for  this  research  was 
9  <  Re  <  110  and  Tong  and  London  had  a  much  larger  range  of 
S  <  Re  <  100,000. 


The  data  for  this  research  did  not  have  the  same 
functional  form  as  Tong  and  London's  data,  thus  a  new 
function  was  sought.  The  new  curve  fits  for  unrolled  screen 

are 

liiy„  =  -9  7  +  3  7to(Re)-05l(to(Re))^  (250  mash).  (4.6) 

la;„  = -II  2-^4  6MRe)-0  63(tn(Re))'  (325  mesh).  (4.7) 

apprrtMingUfjn  Attempts  to  normalize  as  was 

done  with  friction,  /,  t#ere  not  as  successful.  HusMrous 
methods  of  normalisat Ion  were  tried  including  ratios  of 
porosity,  pt  total  heat  transfer  area,  At  minimum  free-flow 
area,  A,t  and  cosribinations  of  these  ratios.  The  least 
scatter  in  the  results  resulted  in  the  normalization 

The  Investigator  defined  jff,  as  the  Colburn 
factor  for  unrolled  screen  and  A^  as  the  total  heat  transfer 
area  for  unrollwd  screen.  Figure  4.17  illustratea  the 
reaulta  of  thia  normalization.  Notice  that  the  five  points 
significantly  above  the  other  points  were  all  for  250  mesh 
screen  This  'jf<-urred  at  the  easimum  flow  rate  of 
approB  isuit  e  1  y  2  9a  10'  ^  kg/s  The  inlet  pressures  were  20 
and  to  ate  The  .  near  legreaston  for  all  the  <iata  was 
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Pigur«  4.17.  Mormallsad  Colburn  Factor  va.  Fraction  Roll,  R 


Tha  Mthodoloqy  to  dataraina  friction  factor,  /,  (Bqs.  4.1 
and  4.2)  appllaa  to  datamlna  ualng  Eq.  4.8.  Tha  Colburn 
factor  for  unrollad  scraan,  can  ba  found  with  Bq.  4.6  or 
4.7. 

Coavactlon  Maat  Tranafar  Coalficiant.  All  tha  work 
with  Colburn  factor,  waa  alaad  at  finding  tha  convaction 
haat  tranafar  coafflclant,  h.  Tha  diaanalonlaaa  taraa  ara 
aatraawly  uaaful  bacauaa  thay  can  ba  appiiad  to  au  Itipla 
altuatlona.  Hoieavar,  tha  aManing  of  tha  data  can  aaally 
bacoaM  obacurad .  For  thla  raaaon  h,  for  tha  apac  i  t  u.  caaa 
of  thia  raaaarch,  la  iiluatratad  in  figa  4.18  and  4.19. 


Although  tha  aparifu  nuadrara  r>nly  apply  to  tha  raganaratora 
of  thla  raaaarch,  tha  trand  froa  tolling  tha  at  raan  la 


clear.  At  low  flow  rates  the  convection  heat  transfer 
coefficient  decreased  an  order  of  magnitude  as  R  was 
increased  from  0  to  0.50.  As  the  flow  rata  increased  the 
difference  bet%feen  rolled  and  unrolled  screen  decreased.  In 
fact,  the  results  for  R  -  0.50  showed  h  actually  becoming 
higher.  Therefore,  at  low  flow  rates  rolling  the  screen  was 
detriSMntal  to  h  and  at  higher  flow  rates  the  differences 
became  less  significant.  These  results  are  generalised  in 
the  heat  transfer  discussion. 

Tests .  Tests  «fere  perfonsed  on  an  es4>ty 
tube  that  was  constructed  in  the  same  Mnner  as  the  other 
regenerators  except  that  screen  disks  were  not  added.  The 
tests  «fere  performed  at  the  minimum  and  maximum  flow  rates 
of  this  research  (0.28  x  and  2.9  x  10'^  kg/s, 

respectively).  For  each  flow  rate  the  inlet  pressure  was  10 
and  30  atm.  Thus,  a  total  of  four  tests  were  performed  at 
the  extreme  conditions  of  the  testing.  At  each  fiow  rate 
the  transient  temperature  response  was  recorded  with  the 
saaw  sMithod  as  all  of  the  other  tests. 

The  results  were  independent  ut  inlet  pressure  so  only 
one  test  from  each  flow  rate  is  shown  in  Fig.  4  20.  Figure 
4.20  provides  a  coe^iarison  between  the  empty  tube  s 
transient  response  and  unrolled  250  mesh  response  The  flow 
rate  and  inlet  pressure  were  the  same.  The  tumparison  wes 
done  with  the  unreal  led  250  swsh  regenerstoi  bei  ause  tta 
screens  had  the  smallest  mass  end  total  heat  transtei  area 
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Pigur*  4.18  Convection  Heet  Transfer  Coefficient,  k,  vs. 
Plow  Rate,  IP,  for  250  Mesh.  Uncertainty  ranged  froai  15 
percent  at  the  Icweat  IP  and  increased  to  50  percent  at  the 

highest  IP  for  each  curve . 
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Pigura  4.19  Convection  Heat  Transtei  Coef t k  lent ,  h,  vm 
riow  Rate,  IP,  for  J25  Mesh  Uncertainty  ranged  froe  15 
percent  at  the  lowest  IP  and  increaatKl  tu  SO  percent  at  the 

highest  *P  for  each  <  urve 
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For  this  roasop,  it  would  be  most  likely  to  be  affected  by 
the  regenerator  caae.  Figure  4.20a  ahows  the  lo%rest  flow 
rate.  The  eaq>ty  tube  had  an  outlet  temperature  7,  which 
remained  roughly  0.5  K  higher  than  the  inlet  temperature, 

T,.  But,  the  important  point  ia  to  compare  the  alope  of  T, 
va.  time  bettfeen  7  and  8  aeconda  for  the  empty  tube  and 
regenerator.  During  thia  time  interval  the  alope  of  the 
regenerator  (right  hand  graph)  waa  euch  ateeper  than  for  the 
ei4>ty  tube.  The  maximum  shpt  uaed  to  calculate  the  nuBd>er  of 
heat  tranafer  unita,  N^,  for  the  regenerator  waa  meaaured 
within  thia  time  interval.  Thua,  the  regenerator  caae 
cauaed  a  minimal  effect  on  the  alope  amaaureaMnt  for  the 
regenerator.  A  aimilar  compariaon  ia  true  for  the  high  flow 
rate  teata  aho«m  in  Fig.  4.20a.  The  only  difference  in  the 
figure  ia  the  that  the  time  elapaed  after  the  drop  in  T,  ia 
the  important  interval  for  the  alope  ccuapariaon.  Thua,  the 
reauita  of  the  eag>ty  tube  teata  ahowed  that  the  regenerator 
caae  had  a  einimal  effect  on  the  awaaureaMint  of  maximum xiope , 
i.e.  heat  tranafer. 

Hmat  Tranafer  DiacuBaion.  Aa  waa  diacuaaed  earlier, 
the  Colburn  factor,  /„,  and  Reynoida  nuad>er,  Rc,  both 
decreaaed  aa  the  acreeri  waa  rolled.  A  rough  approx  lea  t  Ion 
of  can  be  eade  ualng  a  normalisation  of  Colburn  factor 
( Iq .  4.8).  Unfortunately,  the  differences  between  the 
previous  reauita  of  Tong  and  London  (19‘S7)  did  not  allow  the 
investigator  to  use  theii  method  to  provide  an  improved 
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Flgur*  4.20a.  Conparlson  of  Transient  Response  of  the  Empty 
Tube  and  Unrolled  250  Mesh  Regenerator  at  0.28  x  lO''^  kg/s 
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Figure  4.20b.  Cueparisun  uf  Transient  Response  of  tf\e  Empty 
Tube  and  Unrolled  250  Mesh  Regenerator  at  2.9  x  10~^  kg/s 
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curve  fit.  Furthermore,  an  intuitive  understanding  of  the 
actual  change  in  heat  transfer  abilities  can  be  lost  in  the 
dimensionless  value 

The  actual  convection  heat  transfer  coefficient,  h,  for 
the  test  results  provided  some  insight  to  the  heat  transfer. 
The  observations  from  the  specific  values  ot  h  va.  fV  in 
Figs .  4.18  and  4 . 19  can  be  applied  to  general  cases .  For 
Re^  <  40,  rolled  screen  demonstrated  significantly  reduced 
heat  transfer.  But  for  40  <  Re^  <  100,  the  difference  was 
less  pronounced.  The  trend  in  the  range  of  Reynolds  numbers 
tested  suggested  that  the  heat  transfer  could  actually  show 
negligible  differences  for  rolled  screen  as  Re.,  is  increased 
beyond  100. 

Discussion  Qi  rrigtion  And  Hsail  TtAniiar 

The  effects  of  rolled  screen  on  friction  and  heat 
transfer  were  discussed  separately  in  prior  sections.  Now 
the  combined  effects  of  both  will  be  discussed  together. 
Recall  the  term  Re„  is  the  Reynolds  number  for  internal  flow 
based  on  the  properties  of  unrolled  screen. 

Friction  and  heat  transfer  provided  somewhat  opposing 
conclusions.  First,  the  drag  coefficjsnt  per  screen,  i),, 
remained  mostly  unchanged  for  screens  rolled  to  H  s  0.30. 
Though,  between  R  *  0.30  and  H  -  0.50  an  increase  in  (',,  was 
noticed.  However,  heat  transfer  provided  a  very  different 
conclusion.  Any  asKJunt  of  rolling  caused  noticeable 
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reductions  in  heat  transfer,  especially  for  Re„  <  40.  As  Re^ 
increased,  the  relative  reduction  in  heat  transfer  for 
rolled  screen  became  less  significant.  The  most  interesting 
result  was  for  R  =  0.50  at  the  highest  mass  flow  rates.  At 
these  conditions,  h  was  the  same  as  the  unrolled  screen 
(this  favorable  observation  could  be  negated  by  the 
increased  pressure  loses  at  /{  =  0.50).  In  addition,  for 
40  <  Rfi^  <  100  the  relative  difference  in  the  heat  transfer 
was  smaller  than  the  uncertainty  range  of  the  experiments. 
So,  the  data  curves  show  a  difference  in  heat  transfer  due 
to  rolling.  However,  the  uncertainty  is  large  enough  that 
no  conclusions  can  be  drawn  for  40  <  <  100. 

For  RCo  <  40  heat  transfer  ability  decreased  as  the 
screen  was  rolled.  The  drag  coefficient  per  screen,  Cp,  was 
unchanged  for  R  s  0.30.  So,  to  obtain  a  rcKpiired  heat 
transfer  rate,  more  layers  of  screen  may  be  needed 
(increasing  total  area,  i4).  This  in  turn  would  result  in  an 
increase  in  pressure  drop  across  the  regenerator.  Thus,  in 
the  region  of  low  Re^  and  low  R,  there  appears  to  be  no 
advantage  to  rolled  screens. 

The  relative  loss  in  heat  transfer  ability  due  to 
rolling  diminished  as  Rc„  increased.  For  Rc„  •  100  the  heat 
transfer  becasw  the  sasM  for  roiled  screen  as  for  unrolled 
screen.  In  addition,  the  friction  results  did  not  indicate 
an  increase  in  friction  with  increased  rolling  for  up  to 
0.30  fraction  roll.  Therefore,  tests  in  the  region  Rc.  •  100 
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and  R  i.  0.30  showed  the  friction  and  heat  transfer 
properties  were  unaffected  by  rolling.  Thus  the  data  in  the 
region  Rc^  •  100  and  R  i.  0.30  suggest  that  tests  for 
Re„  >  100  would  be  beneficial. 

The  coodbined  effect  of  heat  transfer  and  friction  when 
considered  with  the  dead  space  of  a  complete  cryocooler  is 
%rorth  evaluating.  This  could  most  likely  be  accomplished 
with  the  aid  of  the  numerical  model  "Regen  3.1“  developed 
by  Gary  and  Radebaugh  (1991)  which  was  mentioned  in  Chapter 
I. 
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This  chsptsr  contains  two  s*ct  i  ons  First,  t  h®  tx^t  t  c>*  i  i  n«» 

of  ths  results  is  ststwd  in  the  conclusion.  Second, 
revcosnendst ion  for  future  research  and  laprovesmnts  ti<  the 
tast  apparatus  are  listed. 

Concluaion 

As  the  screen  was  rolled,  Rc  and  f  decreased  The 
graphs  of  r^,  vs.  showed  that  for  fraction  roll  up  to 
0.30  the  drag  coefficient  was  unchanged.  Thus,  the  friction 
results  suggested  that  the  screen  thickness  can  be  reduced 
by  as  euch  as  30  percent  without  increasing  the  pressure 
drop  for  a  fixed  suiss  flow  rate  and  fixed  nuad>er  of  screens. 

In  the  region  Rc„  <  40,  the  heat  transfer  and  A)  was 
decreased  an  order  of  sMgnitude  as  the  roll  fraction  was 
increased  to  0.50.  In  the  region  40  <  Rc^  <  100,  the 
reduction  in  and  h  due  to  rolling  becaac.  less  pronounced. 

Cosd>ining  the  trends  of  the  heat  transfer  and  friction 
results  suggests  that  for  Re„  >  100  rolled  screen  stay 
per fore  the  sasw  as  unrolled  screen.  But,  tests  were  not 
perfonsed  for  Re„  >  100  due  to  lieitations  in  the  present 
test  apparatus.  Thus,  future  research  for  Rt„  >  100  could 
prove  rolled  screen  to  be  a  practical  approach  to  reducing 
the  void  volusM. 
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R*<‘ 


mcl«d  ar«»a^  f  ’r  ’ 


P«rtor»  teatu 


diakMtci  r«<|«rverat  ora  Maintair-  r  h«>  aaai^^  !  lv'%  ra*. 
■o  that  aost  of  the  saaie  apparatus  say  he  used  and  'he 
halxua  usa9e  My  be  ainLeiied. 

2.  Apply  the  friction  and  heat  transfer  results  (  rt>« 
this  research  to  a  Stirling  cycle  SMsdel  .  This  will  answer 
the  question  of  the  conflicting  gain  from  reduced  void 
voluM  and  loss  frosi  reduced  heat  transfer.  A  possible 
SKKlel  is  'Regen  3.1'  (Gary  and  Radebaugh  1991). 


RecooBsended  improveinents  to  the  test  apparatus  include: 

1.  Improve  the  accuracy  of  the  transient  temperature 
measurements,  since  these  measurements  were  the  limiting 
factor  to  the  heat  transfer  calculations.  This  can  be  done 
with  a  high  speed  temperature  data  acquisition  hardware  and 
software  package. 

2.  More  closely  approximate  a  step  temperature  change 
by  removing  the  thermal  effect  of  the  tubing  in  critical 
areas . 

3.  Increase  the  size  of  the  step  temperature  change  by 
improving  the  ice  bath,  because  the  change  was  only  half  ol 
the  desired  20  K  drop. 
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ABMndix  Al  £r9P^gti9g  Ot  HsUam  AD^  1M  stainless  steel 
Table  A. 1  Properties  of  Helium 


PropartiM  of  H*4  at  P  - 101.325  kPa 


M  -  4.003 


r.  K 

p. 

H,  Pi 

'  s 

►.  m'/t 

kj/(k6  •  K) 

*.  W/(m  ■  K) 

Pr 

100 

0  4878 

97  80 

-07 

20  05 

-06 

5.194 

75M 

-09 

0  690 

150 

0.5252 

12.50 

-06 

58.44 

-06 

5.194 

96.90 

-09 

0670 

200 

0  2459 

15  10 

-06 

61.91 

-06 

5.195 

ll.M 

-02 

0M5 

250 

0  1951 

17M 

-06 

90.20 

-06 

5.195 

19.70 

-02 

0.667 

500 

0  16M 

19.90 

-06 

12.24 

-05 

5.195 

15  50 

-02 

0.667 

550 

0.1594 

22.20 

-06 

15.95 

-05 

5.195 

17.20 

-02 

0.670 

400 

0.1220 

24.50 

-06 

19.95 

-05 

5.195 

18  90 

-02 

0  668 

450 

0.1084 

26.40 

-06 

24.55 

-05 

5.195 

20.50 

-02 

0.669 

500 

0.0976 

28.40 

-06 

29.11 

-05 

5.199 

22.10 

-02 

0.667 

MO 

0.0815 

52.20 

-06 

59.61 

-05 

5.195 

25.10 

-02 

0.666 

700 

0.0697 

55.90 

-06 

51.52 

-05 

5.195 

28.00 

-02 

0.666 

MO 

0.0610 

59.40 

-06 

64.62 

-05 

5.195 

90.70 

-02 

0  666 

900 

0.0542 

42.M 

-06 

78.97 

-05 

5.199 

99.40 

-02 

0.665 

1.000 

0.0488 

46.20 

-06 

94.71 

-05 

5.199 

56.00 

-02 

0.666 

1.100 

0.0445 

49.40 

-06 

11.14 

-04 

5.195 

58.50 

-02 

0.666 

1.200 

0.0407 

52.M 

-06 

12.91 

-04 

5.195 

41.00 

-02 

0.665 

1.500 

0.0575 

55.M 

-06 

14.82 

-04 

5. 195 

49.40 

-02 

0.665 

1.400 

0.0948 

M.60 

-06 

16.62 

-04 

5.199 

45.70 

-02 

0.666 

(Kays  and  London,  1984x284) 


Properties  of  304  Stainless  Steel  at  300K 

р,  =  7900  kg/m3 

с,  »  477  J/kg-K 
k,  »  14.9  W/mK 

(Incropera  and  DeWitt,  1990: AS) 
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Appendix  B;  ^  vs .  maximum  slope 


T^d}le  B .  1  WB.  maximum  slope 


1 

0.929 

130 

3.226 

245 

4.422 

15 

1.121 

135 

3.287 

250 

4.467 

1.286 

140 

3.347 

255 

4.511 

25 

1.432 

145 

3.406 

260 

4.555 

30 

1.565 

150 

3.464 

265 

4.599 

35 

1.687 

155 

3.521 

270 

4.642 

40 

1.801 

160 

3.577 

275 

4.684 

45 

1.908 

165 

3.632 

280 

4.727 

50 

2.010 

170 

3.686 

285 

4.769 

55 

2.107 

175 

3.740 

290 

4.810 

60 

2.199 

180 

3.793 

295 

4.851 

65 

2.288 

185 

3.845 

300 

4.892 

70 

2.373 

190 

3.896 

305 

4.933 

75 

2.455 

195 

3.947 

310 

4.973 

80 

2.535 

200 

3.997 

315 

5.013 

85 

2.612 

205 

4.046 

320 

5.052 

90 

2.687 

210 

4.095 

325 

5.091 

95 

2.760 

215 

4.144 

330 

5.130 

100 

2.832 

220 

4.191 

335 

5.169 

105 

2.901 

225 

4.239 

340 

5.207 

110 

2.969 

230 

4.285 

345 

5.245 

115 

3.035 

235 

4.331 

350 

5.283 

120 

125 

3.100 

3.163 

240 

4.377 

355 

5.321 

Curve  fits  from  the  above  table 


maximum  slope _ a _ b 

[0.92,  2.00)  11.77  2.084 

=  a(maximum  slope)^  [2.00,  2.80)  12.19  2.023  (2.12) 

[2.80,  5.32)  12.38  2.008 

Note:  The  QBaslc  progreLJ  listing  on  the  next  page  generated 
this  table.  It  is  based  on  Eq.  2.7. 
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REM  This  program  finds  the  maximum  slope  of  non-dimensional 
temperature 

REM  for  leunda  =0  Note:  calculation  of  Bessel  uses 
logrithmns 

DECLARE  SUB  bessell  (x,  II#) 

OPEN  "blogl.dat"  FOR  OUTPUT  AS  #1 

tau  -  .001:  deltatau  =  .5 

PRINT  #1,  "'Ntu'",  " 'maxslope' " ,  "'tau'" 

FOR  Ntu  *  10  TO  500  STEP  5 
REM  PRINT  #1,  -'Ntu  Ntu 

slope2  =0:  slopel  =  0 
DO  UNTIL  slope2  <  slopel 

X  =  2  *  SQR(Ntu  *  tau) 

REM  PRINT  X 

bessell  x,  II#:  REM  call  the  modified  bessel 
subroutine 

slopel  =  slope2 

slope2  =  Ntu  ^  2  /  SQR(Ntu  *  tau)  *  II#  *  EXP(-(Ntu 

+  tau) ) 

REM  PRINT  #1,  tau,  slope2 
tau  -  tau  +  deltatau 
REM  PRINT  slope2 
LOOP 

PRINT  #1,  Ntu,  slopel,  tau 
PRINT  Ntu,  slopel,  tau 
NEXT  Ntu 

CLOSE  #1 

SUB  bessell  (x,  Inew#) 

REM  hyperbolic  (modified)  bessel  function:  n=l 
REM  using  logrithms 
Ifactl#  =  0:  fact2#  =  0 
k  =  1 

Inew#  =  0:  lold#  =  1 

DO  UNTIL  ABS(Iold#  -  Inew#)  /  lold#  <  .000001 
lold#  -  Inew# 

Ifactl#  =  lfact2#: 
lfact2#  =  lfact2#  +  LOG(k) 

REM  Inew#  =  Inew#  +  x^  (1  +  2*  (k-1))  /2  ''  (1  +  2*  (k 
-  1))  /  factl#  /  fact2# 

Iratio  =(1+2*  (k-1))*  LOG(x)  -(1+2*  (k-1))* 
LOG(2)  -  Ifactl#  -  lfact2# 

Inew#  »  Inew#  +  EXP (Iratio) 

IF  k  =  1  THEN  lold#  =  1 
k  =  k  +  1 
LOOP 
END  SUB 
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SCREEN  PROPERTIES 


regenerator  estimated  using  screen  count  and  the 
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RESULTS  FOR  250  mesh  R  «  0.143 
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RESULTS  FOR  325  mesh  R  -  0.152 
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RESULTS  FOR  325  mesh  R  =  0.303 
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REPEAT  RESULTS  FOR  250  mesh  R  -  0.286  and  EMPTY  TUBE  TESTS 


Appendix  E:  Data  Acquisition  Program 


'*********************************************************** 

120  '  This  is  the  temperature  data  collection  program 
140 

'***********************★«********************************** 
150  SCREEN  0,0, Ot  KEY  OFFtCLStWIDTH  80 
160  CLEAR,  491521 

165  LCXJITB  25,l:COLOR  0,7: PRINT "-PLEASE  WAIT- i COLOR 
7, Ox PRINT"  Loading  DAS-8  I/O  address  and  thermocouple 
loolcup  table  data"  x LOCATE  1,1 
170  DEP  SE6  »  0 

180  SG  256  *  PBBK(6HS11)  -t-  PEEK(6H510) 

190  SG  -  SG  +  491521/16 
200  DEP  SBG  -  SG 
210  BLQAD  "DAS8.BIN",  0 
213  CLS 

220  OPEN  "DAS8.ADR"  FOR  INPUT  AS  #1 

230  INPUT  #1,  BASADR%  'initialize  fc  declare  CALL 

parameters 

240  CLOSE  #1 

250  DAS8  «  0 

260  FLAG%  «  0 

270  1ID%  «  0  'Mode  0  »  initialization 

280  CALL  DAS8  (MD%,  BASADR%,  FLAG%) 

290  IF  FLAG%  <>0  THEN  PRINT" INSTALLATION  ERROR" 

292  'initialize  variables 

293  ZNAZ%»500xZM%-ZMAX%xDIM  D%(2,ZM%) ,TINE(ZM%) 

295  'Get  gain  setting  of  EXP-16 

300  AV>1000x 'CLS X INPUT  "EXP-16  Gain  setting  (100,200,1000 

etc.)*  "/AV 

301  ' - control  menu - 

302  A6E$(0)-"old"xAGE$(l)«"new" 

305  KET(l)  ON 

306  KET(10}  ON 

310  ON  KEY(l)  GOSUB  2000 
315  ON  KEY(IO)  GOSUB  3000 
320  IF  JU  -1  THEN  CLSxJU-0 
322  LOCATE  5,15 

325  PRINT" filename!  " ;N$,AGE$(WARN) 

330  PRINT! PRINT"T( 1)  »  ";T(0) 

335  PRINT! PRINT "T( 2)  -  ";T(1) 

340  PRINTxPRINT"T(3)  -  "?T(2) 

342  PRINT  "CJC  -  -;CJC+273.15 

345  PRINT! PRINT! PRINT"  FI  -  enter  new  filename" 

350  PRINT"  FIO  -  collect  data  " 

355  ZM%-1 

360  GOSUB  600  '  cold  juntion  and  data  collection 
363  Z%-1 
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365  GOSUB  820  '  convert  data  to  temperature 
370  GOTO  320 

549  'I  reduced  the  sized  of  D%  and  T  from  15  to  2  for  this 
program 

550  ' -  STEP  2:  Initialize  an  integer  array  D%(15)  to 

receive  data  — - 

560  DIM  D%(2)  '16  elements,  one  for  each  EXP- 16  channel 

570  'Also  initialize  a  corresponding  real  array  to  receive 
temperature  data 
580  DIM  T(2) 

590  ' 

600  ' -  STEP  3t  Get  cold  junction  compensation 

temperature  - 

610  'Output  of  CJC  channel  is  scaled  at  24 .4mV/deg.C.  This 
corresponds  to 

620  '0.1  deg. C. /bit.  Dividing  output  in  bits  by  10  yields 
degrees  C. 

630  ' 

640  'Lock  DAS-8  to  channel  #7  (CJC  channel  selected)  using 
mode  1 

650  MD%»1  :  LT%(0)=7  *  LT%(1)  =  7 
660  CALL  DAS8  (MD%,  LT%(0),  FLAG%) 

670  IP  PLAG%  <>  0  'THEN  PRINT  "ERROR  IN  SETTING  CJC  CHANNEL" 
t  END 

680  'Next  get  CJC  data  from  this  channel  using  Mode  4 

690  MD%  «  4  t  CJ%«  0 

700  CALL  DAS8  (MD%,  CJ%,  FLAG%) 

710  'Change  output  in  bits  to  real  temperature 
720  CJC  -  CJ%/10 

730  ' -  STEP  4t  Get  the  thermocouple  data  - 


740  CH%  -  0 

750  FOR  Z%  -  1  TO  ZM% 

755  TIME(Zt) -TIMER 
760  GOSUB  1000 
765  NEXT  Z% 

770  RETURN 

790  '  CHI  -  specifies  DAS-8  channel  that  EXP-16  is 
connected  to  (0-7). 

800  '  D%(15)  -  integer  data  array  to  receive  data  from 

channels . 

810  ' 

820  ' -  STEP  5:  Convert  data  to  volts  and  calculate 

temperature  - 

830  'AV  -  Gain  setting  on  Dipswitch  of  EXP-16  (change  to 
suit) . 

835  CJC  -  CJI/IO 

840  FOR  I  -  0  TO  2 

850  V  »  (D%(I,Z%)*5)/(AV*2048) 

865  VJ»1000*V+.992+(CJC-25)*. 040667  'rem  jeff's  curve  fit 
867  T(I)-273. 15+(-. 0142857+25. 8333*VJ-. 595235*VJ^2) '  curve 
fit  temperature 


880  NEXT  I 
890  RETURN 

1000  ' -  Subroutine  to  convert  EXP-16  channels  to  number 

of  bits  - 

1010  'First  lock  DAS-8  on  the  one  channel  that  EXP-16  is 
connected  to. 

1020  LT%(0)  =  CH%  :  LT%(1)  =  CH%  ;  MD%  =  1 
1030  CALL  DAS8  (MD%,  LT%(0),  FLAG%) 

1040  IF  FLAG%  <>  0  THEN  PRINT  "ERROR  IN  SETTING  CHANNEL"  : 
END 

1050  'Next  select  each  EXP- 16  channel  in  turn  and  convert 
it. 

1060  'Digital  outputs  OPl-4  drive  the  EXP- 16  sub-multiplexer 
address,  so  use 

1070  'mode  14  to  set  up  the  sub-multiplexer  channel. 

1080  FOR  SUB%  =  0  TO  2  'note  use  of  integer  index  SUB% 

1090  MD%  =  14 

1100  CALL  DAS8  (MD%,  SUB%,  FLAG%)  'address  set 

1110  IF  FLAG%  <>  0  THEN  PRINT  "ERROR  IN  EXP-16  CHANNEL 
NUMBER"  :  END 

1140  MD%  =4  'do  1  A/D  conversion 
1150  CALL  DAS8  (MD%,  D%(SUB%,Z%),  FLAG%) 

1160  IF  FLAG%  <>  0  THEN  PRINT  "ERROR  IN  PERFORMING  A/D 
CONVERSION" 

1180  NEXT  SUB% 

1190  'All  done  -  return  fron  subroutine 
1200  RETURN 
1210  ' 

2000  ' -  subroutine  to  enter  filename 

2002  IF  WARN  =  0  THEN  2010 

2005  INPUT  "are  you  sure?  (y/n)";A$ 

2007  IF  A$<>"Y"  AND  A$<>"y"  THEN  2120 

2008  CLOSE  #1: CLOSE  #2: CLOSE  #3 
2010  LOCATE  18,10 

2015  INPDT"enter  ' filename' .dat" ;N$ 

2020  NA$="a: "+N$+" .dat" 

2030  NB$="b: "+N$+" .prn" 

2040  NC$="\thermo\data\"+N$+" .dat" 

2050  OPEN  NC$  FOR  APPEND  AS  #1 
2060  OPEN  NB$  FOR  APPEND  AS  #2 
2070  OPEN  NA$  FOR  APPEND  AS  #3 
2080  PRINT  #1, "'filename  =  ";NC$ 

2090  PRINT  #2, " 'filename  =  ";NB$ 

2100  PRINT  #3, " 'filename  =  ";NA$ 

2110  WARN  =  1 
2115  JU=1 
2120  RETURN 

3000  '  -  subroutine  for  data  collection 

3005  LOCATE  18,10 

3010  IF  WARN  =  0  THEN  JU=1:CLS: PRINT  "You  forgot  to  enter 
new  filename.”:  GOTO  3140 
3015  KEY(l)  STOP:KEY(10)  STOP 
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3020  PRINT  "collecting  data...  " 

3030  ZM%=ZMAX% 

3040  GOSUB  600  '  collect  data 
3042  PRINT  "saving  data  ..." 

3044  PRINT  #1,  "'z%'", "'delta 

t' ", " 'Temp(l) ' " 'Temp(2) ' " 'Temp(3) ' " 

3045  FOR  Z%  =  1  TO  ZM% 

3047  GOSUB  820  '  calculate  temperatures 

3048  PRINT  #1,  Z%,TIME(Z%)-TIME(1),T(0),T(1),T(2) 
3050  NEXT  Z% 

3052  CLOSE  #1 

3054  PRINT  #2,  " 'z%' "," 'delta 

t' ", " 'Temp(l) ' " 'Temp(2) ' " 'Temp(3) ' " 

3055  FOR  Z%  =  1  TO  ZM% 

3057  GOSUB  820  '  calculate  temperatures 

3058  PRINT  #2,  Z%,TIME(Z%)-TIME(1) ,T(0) ,T( 1) ,T(2) 
3060  NEXT  Z% 

3062  CLOSE  #2 

3064  PRINT  #3,  "' z% 'delta 

t' ", " 'Temp(l) ' " 'Temp(2) ' ", "'Temp(3) ' " 

3065  FOR  Z%  =  1  TO  ZM% 

3067  GOSUB  82C  '  calculate  temperatures 

3068  PRINT  #3,  Z%,TIME(Z%)-TIME(1) ,T(0) ,T( 1) ,T(2) 
3070  NEXT  Z% 

3072  CLOSE  #3 
3120  WARN  =  0 

3130  KEY(l)  ON:  KEY  (10)  ON 

3135  JU=1 

3140  CLOSE  #1 

3150  CLOSE  #2 

3160  CLOSE  #3 

3180  'ERASE  0%,TINE 

3190  'DIM  D%(2,ZM%) ,TIME(ZM%+1) 

3195  SUB%=0:I=0:ZM%=1:Z%=1 
3200  RETURN 
60000  END 
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Mesh  Size:  230  325  Percent  Roll;  0  15  30  50 

Tank  Pressure  (psig) _  Room  Temp  (  ) _ 

Barometric  Pressure  (mm  Hg) _ 

Target  inlet  Press  (atm)  10  20  30 

Target  Voltage  =  { 14.6962*(atm)-mm  Hg*(0.019337)}/49  965  = _ 

Actual  Voltage _  Actual  Pressure  (Pj)  = _ 

Differential  Pressure 

Actual  Voltage _  AcOud  AP  = _ 

Initial  Temperatures  (degrees  C) 

Ti= _  T2  = _  T3  = _ 

Flowmeter 

Reading _  PmeterO'oi^)” _  Pmeter(Pa)“ _ 

(kg/sec)  = _  Comments: 


Calculated  Results 

(^asic  FUe: _  Grapher  File:  _ Quattro  Pro  Ffle; 


Gtaphermax 


conversion  actor - 


QPro  max—  —  ’ 
dff^aecj 

max  non-dim  slope  = _ 


Table  Ntu  > _ _ - _ )  - 

A(in2)  - _ Ac(m2)  - _ D(m)  =• _ L(m)  = _ p  =• _ 

Re  -  StPr"  =  ^0.667^  = 

-  .  ^  - 
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